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ABSTRACT 

We investigate the launching of jets and outflows from magnetically diffusive accretion disks. Using 
the PLUTO code we solve the time-dependent resistive MHD equations taking into account the disk 
and jet evolution simultaneously. The main question we address is which kind of disks do launch jets 
and which kind of disks do not? In particular, we study how the magnitude and distribution of the 
(turbulent) magnetic diffusivity affect mass loading and jet acceleration. We have applied a turbulent 
magnetic diffusivity based on a-prescription, but have also investigate examples where the scale height 
of diffusivity is larger than that of the disk gas pressure. We further investigate how the ejection 
efficiency is governed by the magnetic field strength. Our simulations last for up to 5000 dynamical 
time scales corresponding to 900 orbital periods of the inner disk. As a general result we observe a 
continuous and robust outflow launched from the inner part of the disk, expanding into a collimated 
jet of super fast magneto-sonic speed. For long time scales the disk internal dynamics changes, as due 
to outflow ejection and disk accretion the disk mass decreases. For magneto-centrifugally driven jets 
we find that for i) less diffusive disks, ii) a stronger magnetic field, iii) a low poloidal diffusivity, or 
a iv) lower numerical diffusivity (resolution), the mass loading of the outflow is increased - resulting 
in more powerful jets with high mass flux. For weak magnetization the (weak) outflow is driven 
by the magnetic pressure gradient. We consider in detail the advection and diffusion of magnetic 
flux within the disk and we find that the disk and outflow magnetization may substantially change 
in time. This may have severe impact on the launching and formation process - an initially highly 
magnetized disk may evolve into a disk of weak magnetization which cannot drive strong outflows. 
We further investigate the jet asymptotic velocity and the jet rotational velocity in respect of the 
different launching scenarios. We find a lower degree of jet collimation than previous studies, most 
probably due to our revised outflow boundary condition. 

Subject headings: accretion, accretion disks - MHD - ISM: jets and outflows - stars: pre-main 
sequence - galaxies: jets - galaxies: active 



1. INTRODUCTION 

Jets as highly collimated beams of high velocity mate- 
rial and outflows of comparatively lower degree of colli- 
mation and lower speed are an ubiquitous phenomenon 
among astrophysical objects. Jets are powerful signs of 
activity and are observed over a wide range of luminos- 
ity and spatial scale. Among the jet sources are young 
stellar objects (YSO), micro-quasars, active galactic nu- 
clei (AGN), and most probably also gamma ray bursts 
dFanaroff fc Riley||1974| | A bell fe Margon||1979] i Mjrnd t] 
& Fried|l98^|Rh'oads|l997||MirabeI fc Rodriguez||l994|). 
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The common models of launching, acceleration, and coili- 
mation work in the fra mework of magnetohydrodynamic 
(MHD) forces (see e.g. |Blandford fc Payne |1982) |Pudritz| 



fc Norman||l983 Uchida~fc~Shibata 1985), although the 



details of the process are not fully understood. 

Jets and outflows from YSO and AGN affect their en- 
vironment, and, thus, the formation process of the ob- 
jects they are launching them. Numerous studies inves- 
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tigate effects of such feedback mechani sms in star for- 
mation and galaxy formation (see e.g. Ba nerjee et al 
2007] ICarroll et aLl[2009| |Gaibler et aT^[20lT[) . However , 
a quantitative investigation of how much ot mass, mo- 
mentum, or energy from the infall is actually recycled 
into a high speed outflow needs to resolve the innermost 
jet-launching region and to model the physical process of 
launching directly. This is the major aim of the present 
paper. 

According to the current understanding, accretion and 
ejection are related to each other. One efficient way to 
remove angular momentum from a disk is to connect it 
to a magnetized outflow. This has been motivated by 
the observed correlation between sig natures of accretion 
and ejection in jet s ources (see e.g 



Cabrit et al. 1990 



Hartigan et al. 1995). 

The overall idea is that the energy and angular mo- 
mentum are extracted from the disk by an efficient mag- 
netic torque relying on a global, i.e. large-scale magnetic 
field threading the disk. If the inclination of the field 
lines is sufficiently small, magneto-centrifugal forces can 
accelerate the matter along the field line. Beyond the 
Alfven point also Lorentz forces contribute to the accel- 
eration. The collimation of the outflow is thought to be 
achieved by magnetic tension due to a toroidal compo- 
nent of the magnetic field. Still, we have to keep in mind 
the fact the the toroidal field pressure gradient acts de- 
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collimating, and an existing external gas over-pressure 
may contribute to jet collimation. 

Before going into further details, we like to make clear 
that with jet formation we denote the process of acceler- 
ating and collimating an already existing slow disk wind 
or stellar wind in to a jet beam. With jet launching we 
denote the process which conveys material from radial 
accretion into a vertical ejection, thereby lifting it from 
the disk plane into the corona, and thus establishing a 
disk wind. 

A vast literature exists on magnetohydrodynamics jet 
modeling. We may distinguish i) between steady-state 
models and time-dependent numerical simulations, or ii) 
between simulations considering the jet formation only 
from a fixed-in-time disk surface and simulations consid- 
ering also the launching process, thus taking into account 
disk and jet evolution together. 

Steady-state modeling have mostly followed the self- 
similar Blandford fe Payne approach (e.g. (Sauty & 



Tsinganos|1994 Contopoulos & Lovelace 1994])), but also 
fully two-dimensional models were proposed ( jPelletier fc 



Pudritz 1992 Li 19931), some of them taking even into ac- 
count the central stellar dipole (Fendt et al. 111995 Paatz 



fe Camenzind|[l996 ). Fur ther, some num erical solutions 



have been proposed by ,(|Konigl et al.| |2010| ISalmeron 

My 



et al. 2011 Wardle & Kohigl 1993 ) m a weakly ionized 
protostellar accretion discs that are threaded by a large- 
scale magnetic field as a wind-driving accretion disk 
They have studied the effects of different regimes for am- 
bipolar diffusion or Hall and Ohm diffusivity dominance 
in these disk. (Self-similar) steady-state mode ls have also 



been applied to the jet launching domain ( Ferreira & 
1995) |Li|1995| [Casse fc Ferreira|2000p , connect 



Pcllcticr 



ing the collimating outflow with the accretion disk struc- 
ture. In addition to the steady-state approaches, the 
magneto-centrifugal jet formation mechanism has been 
subject of a numb er of time-dependent num eric al stud- 



ies. In particula r, Ustyugova et al. (1995) and Ouyed 



& Pudritz (1997) have demonstrated for the first time 



the feasibility of the MHD self-collimation property of 
jets. We note, however, that it was already 1985 when 
the first jet formation simulations were publis hed, in that 



case for much shorter simulation time scale (Shibata & 
Uchida 1985). Among these works, some studies inv es 



tigated artificial collimation ( Ustyugova et al.||1999|), a 
more consistent disk boundary condition ( |Krasnopolsky| 
et al. 1999|), the effect of magne tic diffusivity on colli- 
mation (Fendt & Cemeljic 2002), or the i mpact of th e 



disk magnetization profile on collimation (Fendt 2006). 
In the aforementioned studies, the jet-launching accre- 
tion disk is taken into account as a boundary condition, 
prescribing a certain mass flux or magnetic flux profile in 
the outflow. This is a reasonable setup in order to inves- 
tigate jet formation, i.e. the acceleration and collimation 
process of a jet. However, such simulations cannot tell 
the efficiency of mass loading or angular momentum loss 
from flux from disk to jet, or cannot help answering the 
question which kind of disks do launch jets and under 
which circumstances. 

It is therefore essential to extend the jet formation 
setup and include the launching process in the simula- 
tions. Numerical simulations of the MHD jet la unching 
from accretion disks have been presented first by |Kudoh| 



et al. (1998) and Casse & Keppens (2002), treating the 
ejection of a collimated outflow out of an ey olving-in time 
resistive accretion disk. Zanni et al. ( 2007 1 further devel- 
oped this approach with emphasis on how resistive affect 
the dynamical evolution . An additional cent ral stellar 
wind was considered by Meliani et al. (2006). Further 



studies were concerned about the effects of the absolute 
field strength or the field geometry, in particular inves- 
tigating field strengths around and below equipartition 
i|Kuwabar;i et al.||2005} |Tzeferacos et al.||2009"{ |Murphy 
|et al.||2010[ ). These latter were long-term simulations tor 
several 100s of (inner) disk orbital periods, providing suf- 
ficient time evolution to also reach a (quasi) steady state 
for the fast jet flow. 

Finally, we like to emphasize the fact that jets and out- 
flows are observed as bipolar streams. Jet and counter 
jet appear typically asymmetric with only very few ex- 
ceptions. One exception is the protostellar jet HH212 
showing an almost per fectly symmetric bipolar structure 
(Zinnecker et al. 1998). So far, only very few numerical 
simulations investigating the bipolar launching of disk 
j e ts have been performed. A mo ng them are the works 



of von Rekowski et al. (2003) or von Rekowski & Bran 
denburg| ( |2004[ ) which even included a disk dynamo ac- 
tion. Recent publications consider asymmetric ejections 
of stellar wind com ponents from an offset multi-pole stel- 
lar magnetospher e ( Long et al.|2008 Lovelace et al. 2010 
Long et al. 2012 1, it is therefore interesting to mvesti- 
gate the evolution of both hemispheres of a global jet-disk 
system in order to see whether and how a global asym- 
metry in the large-scale outflow can be governed by the 
disk evolution. This has not been done so far. 

In a series of two papers, we will address both the de- 
tailed physics of MHD jet launching (paper I) and the 
bipolarity aspects of jets (paper II). In the present pa- 
per (paper I) we investigate the details of the launch- 
ing, acceleration, and collimation of MHD jets from re- 
sistive magnetized accretion disks. We investigate how 
the mass and angular momentum fluxes depend on the 
internal disk physics applying long-lasting global MHD 
simulations of the disk-jet system. We hereby investi- 
gate different resistivity profiles of the disk. This paper 
is organized as follows. Section 2 is dedicated to MHD 
equations and to describe the numerical setup, the ini- 
tial and boundary conditions of our simulations. The 
general evolution of jet launching and the physical pro- 
cesses involved are presented in section 3 with the help 
of a reference simulation. Section 4 is then devoted to a 
parameter study comparing jets from different setups. 

In paper II we will present the bipolar jets simulations, 
discussing their symmetry properties and how symmetry 
can be broken by the intrinsic disk evolution. 



2. MODEL SETUP 

We model the launching of a MHD outflow from a 
slightly sub-Keplerian disk, initially in pressure equilib- 
rium with a non-rotating disk-corona. 

As illustrated in Fig. [I] matter is first accreted radially 
along a disk surrounding a central object and is then 
loaded on to the magnetic field lines. The large-scale 
magnetic field is threading the disk and thereby connects 
the accretion and ejection processes. 

Our main goals are 
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Fig. 1. — Cartoon display of the outflow launching process from 
accretion disks. Matter (dashed lines) is accreted along the disk 
surrounding a central object and is loaded on to the magnetic field 
lines (solid lines). The emerging disk wind is further accelerated 
and collimated into a high velocity beam (jet formation). 

i) to determine the relevant mass fluxes (mass ejec- 
tion and accretion rate) , and to study the influence 
on them by the leading physical parameters such as 
magnetic diffusivity, magnetic field strength, 

ii) to determine the resulting jet geometry - that is 
the asymptotic jet radius and opening angle, along 
with the size of the jet launching area of the disk, 
and the asymptotic jet velocity. 

We will in a follow-up paper further extend our setup 
into two hemispheres, to investigate iii) the launching 
of bipolar jets and their symmetry characteristics. The 
relevant mass and energy fluxes for accretion and ejection 
are of essential interest for feedback mechanisms in star 
formation or galaxy formation. We our highly resolved 
simulations of the innermost regions of these objects we 
intend to quantify these properties for a range of possible 
parameters. 

2.1. MHD equations 

For our num erical simulations, we apply the MHD code 
PLUTO 3.01 ( |Mignone et al.||2007| ), solving the conser- 
vative, time-dependent, resistive, mviscous MHD equa- 
tions, namely for the conservation of mass, momentum, 
and energy, 

°" ^ ^ (pv) = 0, (1) 
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Here, p is the mass density, v the velocity, P the thermal 
gas pressure, B the magnetic held, and $ = —GM/R the 
gravitational potential of the central object of mass M, 
withe the spherical radius R = \Jr 2 + z 2 . In general, the 
magnetic diffusivity is defined as a tensor fj (see §2.5). 



The evolution of the magnetic field is described by the 
induction equation, 



dB 
~dt 



— V x (v x B — rj : j) = 0, 



(4) 



with the electric current density j given by Ampere's law 
j = (V x B)/4tt. The cooling term A can be expressed in 

terms of Ohmic heating A = ^r, with T = (fj : j) ■ j, and 
with g measuring the fraction of the magnetic energy that 
is radiated away instead of being dissipated locally. For 
simplicity, here we adopt g = 1. The gas pressure follows 
an equation of state P = (7 — l)u with the polytropic 
index 7 and the internal energy density u. The total 
energy density is 

8tt 



P 



7-1 



pv 



(5) 



Our simulations are performed in axisymmetry apply- 
ing cylindrical coordinates. The CEN03 algorithm as 
3 rd orde r interpolation scheme is used f or spatial inte- 
gration ( Del Zanna & Bucciantini 2002 1 together with 
a 3 rd order Runge-Kutta scheme for time evolution and 
an HLL Riemann solver. For the magnetic field evolu- 
tion we apply the constrained transport method (FCT) 

ensuring solenodality V • B = 0. 

2.2. Units and normalization 

Throughout the paper distances are expressed in units 
of the inner disk radius r- u while pd,i and pd.i denote 
the disk pressure and density at this radius, respectively. 
The index 'i' refers to a number value at the inner disk 
radius at z — and time t = 0. 

In the Appendix \K\ we show for comparison the astro- 
physical scaling for yound stellar (YSO) jets and AGN 
jets. Typically, n < 0.1 AU for YSO and n < 10 
Schwarzschild radii for AGNs. Naturally, we cannot 
treat any relativistic effects of AGN jets with our non- 
relativistic setup. Velocities are measured in units of the 
Keplerian speed ux,i at the inner disk radius. Time is 
measured in units of tj = ri/uK,ij which can be related 
to the Keplerian orbital period tk,i = 27rf;. Pressure is 
given in units of pd.i — £ 2 Pd,iVK i- The magnetic field 
is measured in units of B\ = B z {. As usual we define 
the aspect ratio of the disk e as the ratio of the isother- 
mal sound speed to the Keplerian speed, both evaluated 
at disk mid plane, e = c s /wj<0 For a more details see 
appendix |A"| 



2.3. 



disk 



ana corona 



Initial setup 

We define initial conditions following a setup applied 
by other authors before - a magnetically diffusive accre- 
tion disk is prescribed in sub-Keplerian rotation, above 
which a hydrost atic corona in pressure balance with th e 
disk is located ( |Zanni et aT1|2007| |Murphy et al.||2010[) . 
The coronal density is chosen several orders of magni- 
tude below the disk density, thus implying an entropy 
and a density jump from disk to corona. However, con- 
trary to previous authors (Zarmi et al.||2007 Tzeferacos 
et al. 2009 Murphy et al. 2010 ) we do not apply an initial 



' In PLUTO the magnetic field is normalized considering 4tt = 1. 
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vertical velocity profile in the disk, and only prescribe a 
radial velocity profile. We have seen that for our long- 
term simulations the whole disk system does adjust to 
a new dynamical equilibrium which does not depend on 
the vertical profile of the initial velocity distribution. 

2.3.1. Initial disk structure 

We prescribe an initially geometrically thin disk with 
e = H/r = 0.1 which is itself in vertical equilibrium 
between thermal pressure and gravity. 

We follow the standard setup employed in a number o f 
previous papers ( |Zanni et al.||2007| |Murphy et al.||2010| ). 

As initial disk density distribution we prescribe 



Outflow boundary (Standard PLUTO) 



Pd = Pd,i 
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(Murphy et al. 2010), 



while the initial disk pressure distribution follows 
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The disk is set into slightly sub-Keplerian rotation ac- 
counting for the radial gas pressure gradient and advec- 
tion, 
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following Murphy et al. (2010), but neglecting viscosity. 

In our setup the initial poloidal velocity is imposed by 
hand following the prescription of Zanni et al. (2007), 
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Simulations with initially vanishing disk accretion re- 
sulted in the same asymptotic inflow-outflow evolution. 
Starting with a disk inflow as in Eq. [9j the inflow-outflow 
structure is established on a shorter time scale. 

2.3.2. Structure of the coronal region 

Above the disk, we define an initially hydrostatic den- 
sity and pressure stratification (the so-called "corona"), 
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The parameter 6 = p c /Pd quantifies the initial density 
contrast between disk and corona. We adopt 5 — 10~ 4 . 

2.3.3. Magnetic field distribution 

The initial magnetic fie ld is prescribed by t he magnetic 
flux function tp following Zanni et al. ( 2007 1 , 
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Here B z q measures the vertical field strength at (r = 
n, z = 0). The magnetic field components are calculated 
by rB z = d^/dr and rB r = dV/dz. 
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Fig. 2. — Computational domain consisting of two grids and a 
set of different boundary conditions. The computational domain 
covers (1500 X 1200) uniform grid cells on a physical domain of 
(r X z) = (96.0 X 80.0)n- Another 2000 stretched grid cells are 
attached in vertical direction 80.0 < z < 288.0. For the high 
resolution simulations a physical domain of (50. Ox 50.0)n is covered 
with (2000 X 2000) grid cells to which another 2000 stretched grid 
cells are attached in vertical direction between 50.0 < z < 180.0. 

2.4. Numerical grid 

The computational domain spans a rectangular grid 
region applying a purely uniform spacing in radial di- 
rection and a uniform plus stretched spacing vertically 
(Fig. [2]) . We have run simulations in two different reso- 
lutions (see Tab. fl]). For the low resolution simulations 
the domain extends to (96 x 288) inner disk radii r; on 
a grid of (1500 x 3200) cells, resulting in a resolution of 
Ar = 0.064. For the high resolution simulations the reso- 
lution is increased to 0.025 on the cost of being limited to 
a smaller domain of (50 x 180)ri. For all simulations pre- 
sented here, a uniform grid is used for the magnetically 
diffusive disl0 

For the typical disk model applied with e = 0.1 the 
low resolution provides only 2 grid cells per disk scale 
height at the inner disk radius, while for somewhat larger 
radii, say at r = 5, we have about 10 cells per disk scale 
height. In the high resolution runs the disk r esolution 
is increased by a factor 2.5. As discussed by |Murphy| 



et al. (2010) the low resolution barely resolves the jet 
launching disk surface layer of the very inner disk, where 
steep gradients in density, pressure and diffusivity are 
present. Here, numerical diffusivity is supposed to play 
a role and we have devoted a section to investigate this 
effect. 

In order to follow the jet outflow over long distances 
and in order to provide a sufficient mass reservoir for disk 
accretion a large domain would be wishful. In order to 
resolve the wind launching high resolution would 

6 Although previous publications have applied a stretched grid 
also for the part of the d omain enclosing the disk structure (see 
e.g. |Murphy et al.| (|2010| l) , we note that up to the latest version 
of the code, FLU TO 3.1., the treatment of magnetic diffusivity is 
limited to equidistant grids (see the PLUTO user manual) 
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Fig. 3. — Distribution of the turbulent magnetic diffusivity r) p = 
'?P,i/( r ' 2 ) m the disk-jet system for the reference simulation (see 
Eq. |13[ l, Here, the diffusive scale height is H v = 0.4r. 

be required. However, as mentioned before, PLUTO in 
its current version does not allow for diffusive MHD sim- 
ulations in a stretched grid. Furthermore, we experiences 
that for large domains together with stretched grids (im- 
plying an elongated shape of the outer grid cells) the 
code had difficulties to solve the conservative equation 
and the simulations crashed. 

2.5. Boundary conditions 

For the boundary conditions, axisymmetry on rotation 
axis and equatorial symmetry for the disk mid plane are 
imposed. At the upper z boundary, we use the standard 
PLUTO outflow boundary condition (zero gradient), but 
at the outer radial boundary conditi on, we apply a mod- 
ified outflow condition as derived by |Porth et al.| ( |2011[ ) 
in order to avoid artificial collimation. Most essential is 
the internal boundary enclosing the origin which we call 
a sink. 

2.5.1. Internal boundary - a central sink 

We prescribe a sink for the mass flux in the very inner 
region of the domain. The sink is essential for the follow- 
ing reasons. First, the numerically problematic singular- 
ity at the origin can be hidden by this internal boundary. 
Second, the central sink allows to absorb the accretion 
flow through the disk and thus, emulate a central accret- 
ing object. 

The sink is numerically introduced as an internal 
boundary condition, located on the region r < r; and 
z < z s V\ The internal boundary conditions defined at 
the right and top sides of the sink have significant im- 
pact on the simulation. If ill-defined, spurious effects 
arise during the early evolution already. For the right 
side of the sink, we impose accretion conditions, imply- 
ing a zero gradient for pressure, density, and the vertical 
velocity component. For and a first order extrap- 
olation (constant gradient) is imposed which ensures an 
angular momentum decrease across the boundary which 
is essential for accretion. To allow for accretion further, 
and also to minimize feedback from the sink to the do- 
main, we constrain min(0, u r; ) in the ghost cells for the 
radial velocity component. We assume that the magnetic 
flux is not advected into the central object, and we im- 
pose E<f, — 0. The normal component of the magnetic 
field is calculated from the solenodality condition along 
both sides of the sink . 

On top of the sink, we prescribe the initial local value 
for the gas pressure. In order to avoid evacuation of the 

7 z a is defined corresponding to the resolution, at least 4 cells is 
taken as a sink height in any run. 



regions close to the symmetry axis, we impose a density 
of 110% the initial local density. Effectively, this con- 
dition replenishes the mass into the domain near to the 
axis to overcome numerically difficult low densities. 

2.5.2. Outer boundary conditions 

An essential point of our setup is to impose a proper 
outflow boundary at the outer boundary of the domain 
in r-direction avoiding artificial collimation forces. Here 
we have imple men ted a current-free outflow condition 
(see Appendix Bl ) which avoids spurious collimation 



by Lorentz forces and which has been thoroughly tested 
in our previous papers fPorth fc Fendt|2010 Porth et al. 
2011 Vaidya et al.|2011 ). In order to enable a long-term 
disk-jet simulation requiring a long-living disk accretion, 
we need to provide a sufficiently long lasting mass reser- 
voir. This could be realized by either i) a prescribed 
mass inflow, by a ii) local mass replenishment, or by iii) 
providing a large mass reservoir for jet launching part 
of the accretion disk by extending the outer disk radius 
to large radii. All these approaches have been chosen 
in the literature. We decided to follow option iii) and 
provide a sufficiently high mass reservoir outside of the 
inner launching region by extending the computational 
grid (and thus the outer disk radius) up to about 100 
in ner disk radii. Similar approaches haven been applied 
by |Murphy et al.| ( |2010[ ) ; |Long et al.| ( |2012[ ) . 

We will see below (see Sect. 3.2) that our approach 
is working well, but has, however, its limits if applied 
for a very long lasting simulations. During about 5000 
dynamical time steps we loose about 30-40% of disk the 
mass due to ejection and accretion and due to unwanted 
mass loss across the outer disk boundary. 

2.6. The disk magnetic diffusivity 

A dissipative effect is required for steady-state accre- 
tion in order to allow matter to diffuse across through 
the magnetic field threading the disk. Magnetic diffu- 
sivity allows the mass flux to cross the field lines and 
thus to allow for accretion along the disk. We con- 
sider the magnetic diffusivity to be turbulent in nature, 
or "anomalous", and, thus, cannot be calculated sclf- 
consistently in our setup. The origin of the turbulent 
magnetic diffusivity is usually referred to be as caused 
by the magneto-rot ational instability (MPJ) (Balbusfe 



Hawley||l99ll|l998[ ) in a moderately magnetized disH 8 ) 
iNevertheless, we argue that as we load turbulently m 



if- 

fusive disk plasma into the outflow, it is natural to as- 
sume that the outflowing material is initially turbulent 
initially as well. Therefore, we have the option of ap- 
plying a diffusive scale height H v = e^r larger than the 
thermal scale height H — er of the disk, supposing that 
the turbulent disk material lifted into the outflow will 
remain turbulent a few more scale lengths until turbu- 
lence decays. It is clear that the strong magnetization 
in the upper layers of the disk will prevent generation 
of turbulence by e.g. the MRI and will also quench the 
turbulence in the outflowing material. Without going 
into detail we may estimate the time scale for the decay 
of the turbulence potentially existing in the outflowing 

8 We will later show that we do not resolve the MRI with the 
numerical resolution applied in our disk-jet setup (see §4) 
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material as follows. The decay of (helical) MHD tur- 
bulence follows a power law Em ~ t~ n , where Em is 
the magnetic energy and the power law index depends 
on further details but is in the range of n = 1/2. ..2/3 
( Brandenburg fe Nordlund||2011[ ). Simulations indicate 
that MHD turbulence decays as fast as the hydr o dy- 



et al. 
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With t v > H/c s = er/c s 



= r 1/2 
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bulence launched from the disk into the outflow of e.g. 
r^(r = 5) > 3, thus in the range of half an orbital pe- 
riod. This must be compared to the kinematic time scale 
for the jet launching. If we consider the propagation of 
the initially slow disk wind over a few disk pressure scale 



heights 



i(r = 5) ~ H/u 



wind 



0.5/0.05 = 10, we 



find that this time is of the order of a few eddy turnover 
times. Also, we will later see that jet launching is a rapid 
process with an outflow being established already after 
a few dynamical time steps. 

Given the essential role of the magnetic diffusivity 
for wind/jet launching, we decided to investigate the 
launching process for different strength and different 
scale height of diffusivity. We have run simulations for a 
variety of combinations of parameter values, but for sake 
of comparison, in this paper, we show results for diffusive 

scale height values, e v = 0.1, 0.2, 0.3, 0.4 

Formally , we apply an a prescription ( |Shakura fe Sun^l 
yaev|1973| , similar also to previous works (see e.g. |Zanmj 



et al. 1 2007 1 ) . We assume the diffusivity tensor to be di- 
agonal with the non-zero components 



V</></> = Vp > Vtr = ih 



V4>- 



(12) 



Here, we have defined a poloidal magnetic diffusivity 
rj p = r/^^ = rj Pt if(r,z), and the toroidal magnetic dif- 
fusivity 7fa = r\ rr = r) zz = ?70,i/(r, z), respectively with 
a function f(r, z) describing the diffusivity profile. 

As it is known from the literature, an an-isotropic mag- 
netic diffusivity is requir e d to obtain stationary solutions 



jWardle fe Konigl||1993| |Ferreira fc Pelletier|p95 
|reira||l997p . Most of the numerical simulations fol 



Fer- 
owed 



that approach and did find a stationary-state from their 
simulations as well (see e.g. ([Zanni et al. | |2007 Tzefera- 
cos et al. 1120091 iMurphy et al!ptllo|) )■ Mote that simula- 



tions of (Jasse & Keppcns ( 2002 , 2004[) apply an isotropic 



diffusivity and reach a steady state as well. 

We define an anisotropy parameter by the ratio of the 
toroidal to poloidal diffusivity components, \ = V^a/VpA- 
For the diffusivity function / several options can be con- 
sidered. For the simulations in this paper (paper I) we 
apply 



f = v A (r, z = 0)H cxp -2 



(13) 



where r) Pt i and rj^j govern the strength of mag- 
netic diffusivity. The Alfven speed va {t, 2 = 0) = 
B z (r,z — 0)/y / p(r,z — 0) , and the disk thermal scale 
height H(r) — c s (r, z — 0)/Ox(r, z = 0) are both cal- 
culated along the mid-plane. The parameter e v mea- 
sures the scale height of the disk magnetic diffusivity, 

9 The r/M effectively governs the diffusion of the poloidal mag- 
netic field, while the r\ rr and r\ zz govern the diffusion of the toroidal 
magnetic field 



similar to the hydrodynamic disk scale height parameter 
e v = H v /r. For the present paper, we decided not to 
evolve the diffusivity profile in time. We find that since 
both the disk scale height H and the Alfven speed va do 
vary as the disk evolves, the strength of diffusivity may 
vary substantially from the initially prescribed value (up 
to a factor 10). This highly non- linear feedback may 
disturb the progress of our simulation such that an arti- 
ficially high diffusivity may be derived which will affect 
the numerical time stepping. A constant-in-time pre- 
scription of diffusivity simplifies our aim of disentangling 
the governing physical processes involved in jet launch- 
ing. 

In the follow-up paper II presenting truly bipolar jet 
launching we will discuss further models for the magnetic 
diffusivity. 

We again stress the importance to provide the reader 
with the actual number values of the magnetic diffusivity 
applied^} If we consider a rj^i ~ 1 together with a max- 
imum Alfven speed in the disk of about 10~ 2 (located 
at the inner disk radius), this gives a maximum disk dif- 
fusivity of about r](r,z)i ~ 0.01 — 0.1, through all the 
simulations. Figure K3] shows the distribution of the mag- 
netic diffusivity. The diffusivity is mainly concentrated 
in the disk, increases with disk radius (consistent with a 
decreasing temperature or ionization degree with radius) , 
and decreases with height resembling the transition from 
a cool disk to an ideal MHD outflow. 

2.7. Main simulation parameters 

The simulations are governed by a set of non- 
dimensional parameters for the following physical quan- 
tities, 

• the magnetic field strength defined by plasma-/3, 
and its initial geometry defined by the "bending" 
parameter m, 

• the magnetic diffusivity, with the three parameters 
V^>,i 1 ^p.i) e v governing its strength and anisotropy, 
and the diffusivity scale height, 

• the initial density contrast between disk and corona 

• the initial disk thermal scale height parameter e = 
H/r. 

An overview over these parameters are shown in Tab. [l] 
(first half) for the simulations presented in this paper. 
The second part of Tab. [T] shows the main dynamical 
quantities resulting from our simulations and will be dis- 
cussed below. Table [2] show the similar numbers for the 
simulations of different diffusive scale height. 

3. MHD JET LAUNCHING 

Before presenting detailed results of our parameter 
study of jet launching conditions, we will first discuss 
the general physical processes involved considering our 

10 We remark that in previous simulations which evolve r) in 
time, only the initial diffusivity distribution was provided, butvnot 
its temporal evolution, and, thus, the magnetic diffusivity distri- 
bution which is act ually involved in p roduci ng the outflow (see e. g. 
|Zanni et aI] | |2007| ; [Tzeferacos et al.|j2009| l; |Murphy et al.|j2010| ). 
bee also our numerical example below. 
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Fig. 4. — Time evolution if the jet-disk structure for reference simulation case 1. Shown is the evolution of the mass density (color) and 
the poloidal magnetic field (contours of poloidal magnetic flux * for the levels 0.01, 0.03, 0.06, 0.1, 0.15, 0.2, 0.26, 0.35, 0.45, 0.55, 0.65, 
0.75, 0.85, 0.95, 1.1, 1.3, 1.5, 1.7) for the dynamical times t = 0,500,1000,2000,3000. The dark lines indicate the slow magneto-sonic 
(dashed), the Alfven (solid), and the fast magneto-sonic (dotted-dashed) surfaces at t = 3000. The electric current lines are shown for 
t = 500 (green solid lines). 



long-term reference simulation (case 1). This simulation 
will then later be compared to simulations applying i) 
different diffusivity profiles, ii) a higher grid resolution, 
iii) and different magnetic field strength (see Tab.[T]). A 
similar setup will then be used to launch bi-directional 
outflows, without the constraint on hemispheric symme- 
try (see paper II). 

3.1. General evolution 

Our reference simulation, case 1 is carried out up to 
t = 5000 dynamical time, corresponding to 796 inner disk 
orbital periods and is thus one of the longest simulations 
considering the MHD jet launching. This huge time pe- 
riod corresponds, however, only to 3 rotations at a radius 
r = 40 r\ and correspondingly less at the outer grid ra- 
dius r = 96 n . Consequently, while the inner part of the 
disk, and thus the outflow evolving from that part, has 
reached a steady state situation, the outer disk corona is 
still dynamically evolving. We will therefore concentrate 
our discussion mainly to the inner disk areas. 

The time evolution of density together with the mag- 
netic field is shown in Fig. [4] Note that we do not show 
field lines integrated from certain foot point radii, but 
magnetic flux contours. This allows to visualize the diffu- 
sion and advection of the magnetic field as a consequence 
of the disk evolution. 

The magnetic diffusivity distribution determines the 
coupling between the plasma and magnetic field, and, 
thus, affects the mass loading into the outflow. Since the 
high diffusivity in the outer part of the disk, the coupling 
is weaker, and, thus, the mass loading less efficient. As a 
general result we observe a continuous and robust outflow 
launched from the inner part of the disk, expanding into 
a collimated jet. 

Figure [5] shows the poloidal velocity distribution in jet 
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Fig. 5. — Accretion-ejection poloidal velocity map at t = 3000 
for reference run case 1. Shown is v p distribution in logarithmic 
scale overlaid with which arrows of the normalized poloidal velocity 
vectors indicating the direction of the flow, v p / \ v p \. 

and disk and the normalized poloidal velocity vectors 
indicating the direction of the mass fluxes. 

The outflow is accelerated to super fast magneto-sonic 
speed (see the magnetosonic surfaces indicated in Fig[4| 
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Fig. 6. — Time evolution of the accretion disk. Shown is the 
radial mass flux (pv r ) in the disk for reference run (case 1) for the 
dynamical times t = Q, 1000, 3000, 5000. 

A bow shock develops at the interface between the out- 
flow and the surrounding hydrostatic corona. As the 
outflow propagates, the initial corona is swept out of 
the computational domain together with the bow shock. 
While the interaction with the ambient gas plays a role 
within the first evolutionary steps, the long-term evolu- 
tion of the outflow - its acceleration and collimation - is 
purely determined by the force balance within the out- 
flow and the physical conditions of the launching region. 

3.2. Disk structure and disk mass evolution 

As jets are launched from disks, the disk evolution is 
itself an essential part of the jet evolution and must be 
carefully considered. It is expected that the jet mass flux 
would somewhat correspond to the disk accretion rate 
(which would also depend on actual mass of the disk). 

Our simulations show that the disk evolution is com- 
plex, with smooth accretion phases are followed by more 
turbulent stages. The time evolution of the accretion 
stream pv r is shown in Fig. [6j We see that accre- 
tion starts first at small radii and fully establishes after 
t = 1 000. Accretion shocks and vortices may destroy the 
quasi-steady state situation. Close to the outer disk ra- 
dius, mass is lost through the outer boundary. At late 
evolutionary stages we see a complex pattern in the pv r 
distribution indicating small shock and turbulent motion 
within the disk. 

Nevertheless, the negative mass fluxes in Fig. |6j to- 



gether with the negative velocity vectors in Fig. [5] clearly 
indicate that accretion is established throughout the 
whole disk. 

We now discuss the time evolution of the disk mass 
for our reference simulation. In order to measure the 
disk mass, we need to carefully define the control volume 
defining a corresponding "disk surface". Here, we con- 
sider as disk surface the location where the mass density 
at each radius falls bellow 10% of the mid-plane density 
(because of the vortex motion we could not use the nega- 
tive velocity as proxy for accretion). Measuring the mass 
contained in subsequent disk rings dM/dr = 2irrpdz, we 
see that most of the mass is indeed stored within the 
outer disk areas (Fig. [7j top). We may identify two 
essential effects affecting the long-term disk evolution. 
Firstly, the mass reservoir in the outer disk has decreased 
- mainly due to outwards mass loss across the outer disk 
boundary. Secondly, the mass of the inner jet-launching 
disk has also decreased - due to both outflow activity and 
accretion into the sink. 

Comparing the mass content of the disk at the initial 
and final stage, we find that the total disk mass decreases 
from M ~ 520 at t = to M ~ 320 at t = 5000 (in 
code units). Integrating the mass loss by accretion and 
ejection into the jet (see Fig. [7]) over 5000 dynamical 
time steps, we find a total mass loss of M ~ 75. The 
difference, i.e. M ~ 125, is the mass which is lost from 
the outermost disk in vertical and radial direction. 

In other words, a substantial amount of the disk mass 
loss is lost from the very outer part of the disk and does 
not directly influences the jet formation. The disk mass 
which is lost from the inner part is partly lost as disk 
wind/jet, partly accreted into the sink and partly re- 
plenished by accretion from outer radii. 

The time evolution of the disk mass gives a similar pic- 
ture (Fig. [TJ bottom). Integrated over the whole domain, 
the disk loses about 38% of its mass until t — 5000 in the 
reference simulation case 1. If we decrease the integra- 
tion domain (outer radius r = 50) the inner disk looses 
less mass while part of its mass is being accreted from 
outer disk radii. In general this implies that up to time 
t = 5000, there is i) still sufficient mass to support the 
accretion process, but also that ii) the disk mass shows a 
considerable decrease which may have changed the inter- 
nal disk properties. Thus, for an even longer-lasting sim- 
ulation one would have to invent another mass source for 
the disk accretion (e.g. by a physical mass inflow bound- 
ary condition at the outer disk radius properly taking 
into account angular momentum transfer, or a floor den- 
sity distribution in the disk). 

Along with the hydrodynamic disk evolution, the mag- 
netic field distribution evolves as well, subject to the 
competing processes of advection and diffusion. This also 
changes the launching conditions, as the overall profiles 
of plasma-/?, resp. magnetization, are affected. Figure [8] 
displays the time evolution of the magnetic flux surface 

= 0.1 for reference simulation case 1, initially rooted 
at the radius r = 2.0. We see that this flux surface is 
first moving (diffusing) outwards "driven" by magnetic 
pressure gradient and tension, and then, once disk accre- 
tion has established again moves inwards being advected 
with the mass flux. After about t = 3000 both processes 
balance and a quasi steady-state situation in the system 
evolution is reached. Considering simply flux conserva- 
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Fig. 7. — Time evolution of the disk mass. Shown is the radial 
profile of the mass distribution dM/dr (integrated mass per disk 
ring, in code units) at times t = 0,3000,5000 (top), and the time 
evolution of the total disk mass normalized to the initial disk 
mass M(j j for integration radii ri = 50, 96 (bottom). 

tion, we may estimate the change in poloidal magnetic 
field strength and the corresponding change in magnetic 
diffusivity and plasma-/?. 

Estimating an average field strength B p and a cor- 
responding magnetic flux 'J ~ B p r 2 , flux conservation 
tells us that * t =iooo ~ V>t=50- Therefore, £? Pi t=iooo ~ 
10 -B p ,t=50: since this flux surface (e.g. ^ = 0.1) is now 
rooted at a different radius. With the increased poloidal 
magnetic field strength the Alfven speed «a ~ B p in- 
creases, and, thus, also the magnetic diffusivity param- 
eterized as rj ~ va- Thus, w e expect that in previou s 
work (e.g. |Zanni et al] ( |2007[ ) ; [Tieferacos et al] ( ]2009[ )) 
the actual (time-dependent) value of the magnetic diffu- 
sivity may differ substantially from its initial value. We 
believe that this will impact the derived mass fluxes. A 
similar estimate can be made for the plasma-/3 or mag- 
netization. Since /3 ~ B~ 2 , the increase in plasma-/? is 
by a factor 100, which also is expected to affect the jet 
formation severely. We will come back to this point later 
when we compare different parameter runs. 



3.3. Accretion rate and ejection efficiency 

The main goal of this paper is to investigate what frac- 
tion of the mass accretion is loaded into the outflow, and 
how the mass loading is affected by various disk param- 
eters. The derived mass ejection-to-accretion ratio can 
be an important ingredient for studies of AGN or YSO 
feedback, or another similar self-regulating outflow sce- 
nario. 

Efficient accretion is feasible only if angular momen- 
tum is sufficiently removed. Since we do not consider 
a physical viscosity in our simulations, angular momen- 
tum removal from the disk is mostly accomplished by 
the torque of the magnetized disk wind. Therefore, in 



our simulations disk accretion can only work if a disk 
wind has been established. In the following we discuss 
the inflow and outflow rates in our disk-jet simulations, 
concentrating first on reference run (case 1). 

In order to calculate the integrated properties of inflow 
and outflow, such as mass flux or angular momentum 
flux, we need to carefully select the integration domairFj 
The derived fluxes depend strongly on the choice of the 
integration boundaries - thus, on how we distinguish be- 
tween material belonging to accretion or ejection. In re- 
ality there is no disk surface, but a smooth transition 
between accretion and ejection. The initial setup of a 
thin disk with aspect ratio e = 0.1 dynamically evolves 
in time, and so does the disk surface. 

Our control volume to measure the disk accretion and 
ejection rates is defined by an axisymmetric sector en- 
closed by two surfaces perpendicular to the equatorial 
plane located at r = 7*1 and r = r 2 , and two other sur- 
faces being the disk mid plane, and a surface Si which 
is close by and parallel to the initial disk surface (see 
Fig. [9]). We usually adopt r% = 1.0. 

The accretion rate is then calculated as 



M acc (r) = -2nr I pv r dz 



iH 



(14) 



where the parameter a controls the scale height of the 
integration, and H = H(r, t = 0) = er is the initial ther- 
mal scale height of the disk. Considering the evolution 
of the large-scale disk velocity (see Fig. pL we have cho- 
sen a = 1.6 for our reference simulation. Similarly, we 
calculate the ejected mass flux as 



cjccMIs 



■dA, 



(15) 



where the integration is done along the inclined surface 
area S from r\ to r 2 considering the area element dA s . 

The time evolution of the accretion and ejec tion rate 
for the reference simulation is shown in Fig. 10 We have 
calculated the accretion and ejection rates for different 
sizes of the integration domain (thus with different outer 
radius r 2 ). We observe that for small radii the accretion 
rate saturates at an approximately constant value. For 
the outer parts of the disk, however, a longer dynamical 
time is required for the disk to evolve into a new dynam- 
ical steady state. This is visible in the time evolution of 
the accretion rates - for increasing radii r 2 = 3, 5, 10, the 
time when the plateau phase is reached, is increasing to 
t = 300, 800, 1500. Beyond the plateau phase, the mass 
fluxes slightly decrease, most probably due to the over- 
all decrease of the disk mass itself, and the subsequent 
change in the internal disk dynamics. The accretion rates 
at large radii are larger than those at smaller radii - the 
mass difference is ejected as outflow. For the control vol- 
ume with larger r 2 the ejection rate increases, which is 
simply a consequence of the inc reased integration area. 

The last panel in Fig. 10 shows the ejection-to- 



accretion rate Af e j ec /M acc . Again, depending on the size 
of the control volume this fraction changes. Although 
the outflow quickly evolves from the disk surface, and 
extends soon to large radii, several orbital periods are 
required to establish full accretion. At earlier times and 

11 Averaging in time is applied for all the fluxes evolutions. 
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n r2 

Fig. 9. — Control volume to measure the accretion and ejection 
rates in the jet-disk structure. Integration is at/between the radii 
7*1 and T2 along the disk. The integration surface denoted by S is 
inclined and is parallel to the initial disk surface. 

for large radii, the accretion process is not fully estab- 
lished, resulting in a somewhat arbitrarily high ejection- 
to-accretion ratios even above unity. For the inner part 
of the disk within radii < 10, about 60% of the accret- 
ing material is diverted into the outflow for our reference 
simulation. This is a large fraction and similar to simula- 
tions in the literatur e, but definitely more than derived in 
steady state m odels ( Pelletier fc Pudritz|19 92; Fcrreira & 
Pelletier|1995[ ) . For the other cases we have investigated, 
also smaller rates were obtained (see below). 

Applying a ra dially self-similar approxim ation of the 
MHD equations, Ferreira fc Pelletierj ( 1995 ) have intro- 
duced a so-called ejection index 4, 



1-1^ 



(16) 



essentially resulting from local mass conservation (where 
r\ = r\ a nd r e = 7 *2 in o ur notation). With self-similar 
solutions Ferreira (1997) constrained the ejection index 
to 0.04 < 4 < U.U8. In comparison, for our reference run 
we find both larger numerical values and also a wider 
range for the ejection index, 0.1 < 4 < 0.5 . Table [T] 
provides an overview over the ejection indexes we mea- 
sure. We can apply our reference run to different astro- 
physical sources (see[A]), thus for a stellar jet with the 
central mass « IMq and p\ w 10~ 10 , the accretion rate 
is ss lO _6 M yr _1 . This value for AGN with the central 
mass « 10 s M© and p- sa 10~ 12 is about w 7M yr _1 . 

3.4. Launching forces and outflow formation 

In order to investigate the launching, the acceleration 
and the collimation processes of the outflow, we now con- 
sider the forces acting in the disk-outflow system. 

To identify the forces for acceleration and collimation 
explicitly, it is preferable to project them along or per- 
pendicular to a certain flux surface, respectively. By 
comparing these projected force components, we may 
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Fig. 10. — Time-evolution of the mass fluxes. Shown is evolution 
of accretion rate (top), the ejection rate (middle), and the ejection- 
to-accretion ratio (bottom) for the reference simulation case 1, all 
in code units. Colors indicate mass fluxes calculated for increasing 
outer radii of the control volume, T2 = 3, 5, 10, 15, 20, 30. 
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Fig. 11. — Accelerating and collimating forces. Shown are the 
profiles of the parallel (top) and perpendicular (bottom) specific 
forces in logarithmic scale for simulation case 1 at time t = 1000 
along the distance A of the field line rooted at (5,0). Here 
Fq, Flu Fp denote the centrifugal force, gravity, the Lorentz 
and the gas pressure gradient forces ( all in absolute value) , respec- 
tively. The vertical lines present the slow-magnetosonic (dashed 
line), the Alfven (solid), and the fast-magnetosonic (dot-dashed) 
point, respectively. 

disentangle the do minant driving and collimation mech- 
anism. Figure [TT] shows the force components along a 
flux surface (hemline) rooted at (5, 0), and also the crit- 
ical MHD surfaces - the slow-magneto-sonic, the Alfven, 
and the fast-magneto-sonic surface. The forces involved 
in driving the outflow are the centrifugal force Fc, the 
Lorentz force F^, the gas pressure gradient Fp, and 
gravity Fq. Among them, the pressure gradient and 
the centrifugal force are de-collimating while gravity and 
Lorentz force have collimating components. The pres- 
sure gradient, centrifugal force and Lorentz force also 
contribute to accelerate the outflow. In agreement with 
the literature we see that upstream the slow magneto- 
sonic point, the acceleration is mainly by the gas pres- 
sure gradient. For the super-slow outflow the Lorentz 
force and the centrifugal force dominate. For the super- 
fast flow the Lorentz force plays a major role. In case 
of the collimating forces the situation is different. Be- 
fore the Alfven point, the de-collimating centrifugal force 
dominates. Beyond the Alfven point the Lorentz force 
has a comparable contribution, and finally after the fast- 
magneto-sonic point, the Lorentz force dominates and 
collimates the outflow. The Lorentz force behavior is a 
sign of the electric current distribution in the disk/jet 
system. It compresses the disk while is collimating the 
corona. 

The forces involved in launching the outflow are mainly 
the vertical gas pressure gradient which counteracts the 
tension forces of the poloidal field and gravity (see 
Fig. 11). The thermal pressure gradient is always posi- 
tive, and thus supports launching. It increases from the 
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Fig. 12. — Launching forces of the outflow. Shown is the vertical 
profile of the vertical components of gravity, the pressure gradient 
force, and the Lorentz force (top), as well as the vertical profile of 
the net force (bottom) at time t = 1000 at a radius r = 5 for in 
case 1. 

disk mid plane to the surface and then decreases in the 
corona. The tension term is mostly negative, thus com- 
pressing the disk, and does not support launching. The 
magnetic pressure (toroidal component and total) are 
both positive above the disk, and, thus, may support 
the launching process. The toroidal component of the 
Lorentz force i^,L provides the magnetic torque brak- 
ing the disk, and loading the outflow with both angular 
momentum and energy (not shown). In order to brake 
the disk, the torque must be n egative in the disk and 
change sign at the disk surface ( |Ferreira 1997). This is 
confirme d by our simulation. 

Figure [12] show for comparison the three main force 
components (top) and the net force (bottom). We see 
that in the disk the vertical force components almost 
balance, while above the disk a considerable net force 
remains which launches and accelerates the outflow. In 
the disk the (positive) gas pressure is almost balanced by 
the (negative) Lorentz force. However, the gas pressure 
slightly wins and counteracts gravity. Essentially the fig- 
ure shows that launching is a process which happens at 
the disk surface. The outflow material is not lifted from 
the mid plane into the disk wind. It is the disk material 
which is accreting along the disk surface which is loaded 
into the outflow. The net (specific) force which is respon- 
sible for launching and initial acceleration is about 10% 
of the value of the single force components. 

3.5. Jet radius and opening angle 

The radius of the asymptotic jet and its opening an- 
gle can be measured by the observations rather easily. 
Therefore it is interesting to provide a comparison from 
the simulations. Simple energy conservation arguments 
tell us that jets with the observed kinetic energy must 
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TABLE 1 

Input parameters and derived dynamical parameters of our simulation runs following a magnetic diffusivity profile Eq. [T3] 

WITH A SCALE HEIGHT OF THE DISK DIFFUSIVITY e v = 0.4. DISPLAYED ARE THE INPUT PARAMETERS OF GRID RESOLUTION Ar, Az, 
MAGNETIC POLOIDAL AND TOROIDAL DIFFUSIVITY r)p,i AND 77^;, THE DIFFUSIVITY ANISOTROPY \, AND THE PLASMA-BETA PARAMETER /3. 
THE FOLLOWING OUTFLOW PARAMETERS ARE DERIVED FROM OUR NUMERICAL SIMULATIONS: THE ASYMPTOTIC JET RADIUS rj ct , CALCULATED 
AS MASS FLUX WEIGHTED (SEE EQ. \17\ ■ THE LAUNCHING RADIUS T\, THE TY'PICAL SPEED OF THE OUTFLOW, l>j et , THE AVERAGE COLLIMATION 

DEGREE Cj e t (SEE EQ. |l8|, THE MASS ACCRETION RATE M acc , AND MASS EJECTION RATE M c j cc , THE EJECTION INDEX £. THE DIFFUSIVE 
SCALE HEIGHT e v THE KINETIC AND MAGNETIC ANGULAR MOMENTUM LOSSES BY THE OUTFLOW, j kin AND j mag , RESPECTIVELY, AS WELL AS 

THE TOTAL ANGULAR MOMENTUM LOSS Jtot , AND ITS ASYMPTOTIC KINETIC LUMINOSITY L kin = 0.5M z V? ct z . ALL PARAMETERS ARE GIVEN 

IN CODE UNITS. ' 





case 1 


case 2 


case 3 


case 4 


case 5 


case 6 
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case 8 
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case 10 


Ar 


0.064 


0.064 


0.064 


0.064 


0.064 


0.025 


0.025 


0.064 


0.064 


0.064 


A -v 


U.UDD 


U.UDD 


U.UDD 


U.UDD 


U.UDD 


U.Uzo 


n no k 


U.UDD 


U.UDD 


U.UDD 


f_ 
>l 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


Vp,i 


0.03 


0.01 


0.15 


0.09 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


V<t>,i 


0.09 


0.03 


0.09 


0.27 


0.01 


0.09 


0.09 


0.09 


0.09 


0.09 


X 


3 


3 


3/5 


3 


1/3 


3 


3 


3 


3 


3 




10 


10 


10 


10 


10 


5000 


10 


50.0 
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r jot,z=180 


43.8 


73.8 


24.46 


28.57 


33.19 






47.8 


50.0 


34.13 


rjct,z=60 


25.01 


38.41 


15.37 


17.72 


25.02 


21.69 


19.6 


27.54 


24.18 


17.22 


n 


3.8 


8.2 


3.0 


4.2 


3.0 




3.5 


2.8 


0.8 


0.7 


fjct,z=170 


0.4 


0.23 


0.85 


0.56 


0.8 


0.47 


0.42 


0.56 


0.49 


0.27 


fict,z=280 


0.45 


0.3 


0.92 


0.61 


0.9 






0.5 


0.89 


0.22 


Cjet.lQ 


2.37 


0.013 


5.28 


5.1 


2.58 


0.89 


4.6 


1.56 


0.32 


2.08 


M z 


0.0005 


6 x 10~ 6 


0.0002 


0.0002 


0.0007 


5 x 10~ 5 


0.0008 


0.0001 


0.0004 


0.0002 


Cjct,2 


2.90 


0.3 


4.45 


6.7 


2.28 






2.46 


0.42 


1.84 


M z 


0.001 


0.0002 


0.0003 


0.0004 


0.001 






0.0003 


0.0006 


0.0004 




0.015 


0.022 


0.0038 


0.0035 


0.022 


0.001 


0.013 


0.005 






Mejec 


0.008 


0.016 


0.001 


0.001 


0.007 


0.005 


0.009 


0.004 






Mcjcc /Macc 


0.5 


0.7 


0.2 


0.2 


0.31 




0.6 










0.3 


0.5 


0.09 


0.09 


0.1 




0.39 








•/kin 


0.008 


0.005 


0.002 


0.002 


0.027 


0.005 


0.005 


0.008 






^mag 


0.03 


0.01 


0.0025 


0.005 


0.034 


0.0 


0.016 


0.008 


0.0 


0.0 


Jtot 


0.038 


0.015 


0.004 


0.007 


0.06 


0.005 


0.021 


0.017 








1 x 10~ 4 


9 x 10" 6 


1.3 x lO" 4 


7.4 x 10~ 5 


4 x lO" 4 


5.5 x 10" 6 


7.1 x 10~ 5 


3.7 x lO" 5 


1.9 x 10~ 4 


10~ 6 



a) ]The collimation degree is measured from two different areas, fj ot i enclosing < r < 40, < z < 160, and £j ct 2 with the area 
< r < 80, < z < 250. 

b) ]The mass flux values for the two simulations case 9 and 10 with high ft are omitted due to their highly perturbed behavior. 

of the outflow. For reference run case 1 we find a mass 
flux weighted asymptotic jet radius for the high-velocity 



originate from a disk area very close to the central object. 
This is a fact which seem to be hold for all astrophysical 
jet sources. We will later consider the question of how 
these asymptotic properties are related to the conditions 
of the jet launching process. Here we first provide a clear 
numerical definition of these properties and apply them 
to our reference simulation. In order to do a quantitative 
comparison, we need to define a "radius of the outflow". 
We suggest a definition of an axial mass flux weighted jet 
radius, measured at a certain distanc e z = z m from the 
source (see also ( |Porth fc Fendt||2010[ )), 



[J^rpv z dr] ; 
[Jo™ P v z dr ] , 



(17) 



This radius measures the bulk of the mass flux contained 
in the jet at a specific distance z m from the mid plane. 
In order to derive the jet launching area, we now adopt 
that flux surface which follows the bulk of the mass flux, 
thus passing through the point (rj et (z m ), z m ), and trace 
the same flux surface back to the disk surface where it is 
rooted. This foot point radius defines the launching area 



component of 



44 at z m = 180, corresponding to 



4AU for a protostellar scaling applying r\ = 0.1 AU. 
For the jet launching radius we measure n ~ 0.4 AU 
in our reference run. Similarly, the launching area for 
the extended low veloci ty component (LVC) of the DG 
Tauri wind measured by Anderson et al. ( 2003 ) extending 
from ~ 0.3 to 4 AU from the star. They have invented 
a method relying on the MHD energy and angular mo- 
mentum conservation along the jet. By comparing the 
observed kinetic energy and jet rotation they infer the 
necessary disk rotation and launching area in the wind- 
launching region. 

The asymptotic jet opening angle is the other jet char- 
acteristics. One way to measure the jet opening angle 
is by the inclination angle of the flux surface defined by 
the asymptotic jet radius as discussed above. A nother 
measu re of the collimation has been suggested by |Fendt| 
( 2006 1 who assigned an average collimation degree £ of 
the outflow by comparing the vertical and radial mass 
fluxes in the outflow (applying a proper normalization 
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Fig. 13. — Comparison of different parameters runs (see Tab. [Tjlat t = 3000. Shown is the density distribution nd the poloidal magnetic 
field (contours of magnetic flux for levels <f = 0.01, Q.03, O.OfjTu.l, 0.15, Q.2, 0.26, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.1, 1.3, 1.5, 1.7) for 
simulation runs with rather low diffusivity (case 2), rather high poloidal diffusivity with anisotropy \ = 3/5 (case 3), rather high diffusivity 
strength (case 4), rather low toroidal diffusivity with anisotropy \ = 1/3 (case 5), simulations with a rather weak magnetic field with 
plasma-/? = 5000 (case 6), plasma-/? = 100 (case 8), plasma-/? = 250 (case 9), and also a simulation with three times higher resolution (case 
7), all to be compared to our reference simulation case 1, see Fig. [5] 



per surface area of a cylinder of height z = z m and ra- 
dius r = r m ), 



Cjct — 2 



z m [J rm rpv z dr] ZiB 

rm [Jo 5 "" r mPVidz] r „ 



(18) 



where only positive poloidal velocities arc considered. 
Outflow collimation would simply imply that £j c t > 1. 
For our reference run, we find (j Ct = 2.37 implying that 
about two times more mass is propagating along the jet 
that away from the jet axis. For comparison, the flux sur- 
face which encloses the bulk of the jet mass flux, has an 
(half) opening angle of about 10° at (r, z) = (44,180), 
but collimates slightly more further downstream. We 
note that this definition is related also to the concentra- 
tion of mass flux across the jet, and, thus, provides a 



somewhat different information than the opening angle 
of the field lines. For example, a cylindrical jets with 
zero degree opening angle may have a narrow or a broad 
radial density or mass flux profile. With our definition 
the narrow mass flux would be interpreted as more col- 
limated. 

4. COMPARISON OF PARAMETER RUNS 

In this section we compare simulations governed by 
different input parameters whic h rely on a magnetic dif- 
fusivity profile following Eq. [l3]with a scale height of the 
disk diffusivity H v = e„r, grid resolution Ar, Az, maxi- 
mum magnetic poloidal and toroidal diffusivity ry Pi i and 
77^4, the diffusivity anisotropy x, and the plasma-beta 
parameter /? (see Tab. [I]). We understand these param- 
eters as the governing parameters of jet launching from 
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Fig. 14. — Comparison of mass fluxes for the different cases. 
Shown is the time evolution of the accretion ejection rate (top) and 
the ejection rate (bottom) for simulation runs applying a different 
parameter setup (see Tab.Q. 



accretion disks. All simulations have been carried out up 
to t = 3000 (reference run up to t = 5000). We found 
that this time is sufficient to reach a quasi steady state 
of accretion-ejection (roughly after t = 1000). 

We first show a comparison of the global mass density 
distribution at t = 3000 for simulation setups resulting 
from a different strength rj or different anisotropy x 01 
diffusivity, a, di fferent plasma-/?, or a different grid reso- 
lution (Fig. 13). The immediate result is that both the 



disk structure and the dynamical evolution of the out- 
flow, change substantially compared to the reference run 
(case 1). Jet-like outflows have been formed in all cases, 
although in some cases like case 6 with /3 = 5000, or case 
9 with j3 — 250 the outflow appears highly filamentary. 
Denser outflows are observed as in case 7 with higher res- 
olution, or case 2 lower diffusivity, or case 5 with lower 
toroidal diffusivity, while in other cases the outflow is 
more tenuous. Outflow and disk evolution is interrelated 
- a denser outflow thus implies to a geometrically thinner 
accretion disk, as more of the accreting matter has been 
diverted into the outflow. The smoothness of the out- 
flow varies for the different setups. We obtain a much 
more filamentary and perturbed structure for outflows 
for which the (physical or numerical) diffusivity in the 
disk is lower. A similar statement can be made for a low 
magnetic field strength. 

We now consider the main properties of the accretion- 
ejection system in quasi steady state for all simulation 
runs. The two most prominent physical quantities are 
the mass and the angular momentum fluxes. Similar to 
the mass fluxes defined in Eqs. [14] and [15] we integrate 
the ejection torque that is the torque exerted on the disk 
by the outflow applying the same control volume (see 
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10 




Fig. 15. — Comparison of angular momentum fluxes for the dif- 
ferent cases. Shown is the time evolution of the vertical angular 
momentum flux from the diskjCalculatcd for a control volume with 
n = 1.0, r 2 = 10.0 (see Tab.[l}. 



§3.3), 

ikin = / rpv.pijp-dAs, j mag 



-^B p -dA, (19) 



with the kinetic and magnetic angular momen tum flux, 
Jkin and J mag , carried by the outflow. Figure 
the time evolution of the accretion rate (top)' 
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shows 
"and the 

ejection rate (bottom) for the different cases. The accre- 
tion rate is calculated for radius r = 10 for all easeful 
For comparison, we show the correspo ndin g vertical an- 
gular momentum fluxes evolution (Fig. [l5]) . For all cases 
investigated, accretion sets in after several 100s of ro- 
tations and is fully established within t ~ 1000. After 
some initial fluctuations, the accretion rate levels off into 
a steady state, depending on the physical parameters pre- 
scribed in the simulation. 

Depending on the efficiency of the angular momentum 
transfer from the disk, the disk establishes a different 
accretion rate. We find that the simulation runs with 

12 For high plasma-/3 (case 9, 10), the mass fluxes evolution looks 
strongly fluctuating and are therefore omitted. 
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the highest accretion rates, also have the hig hest angular 
momentum flux from the disk (Figs. 14|15 1. These are 
the the runs case 2 with rather low diffusivity and case 
5 with rather low poloidal diffusivity, strongly indicating 
that with a low diffusivity and, thus, a strong coupling 
between field and matter, the magneto-hydrodynamic 
torque of the jet on the disk is most efficient, and subse- 
quently leads to a more efficient disk accretion. Similarly, 
the high plasma-/? (thus with a low field strength, simu- 
lations runs (case 6, 8, 9 and 10) exhibit a weak magnetic 
angular momentum removal and, consequently, are ineffi- 
cient accretors (for a summary see Tab.nl). In summary, 
we find that the extraction of angular momentum from 
the disk by the outflow and accretion are clearly interre- 
lated. 

In the following sections we discuss how the govern- 
ing system parameters such as diffusivity, plasma-/?, and 
numerical resolution affect the mass flux evolution. 



4.1. Possible impact of numerical diffusivity 

Before we further investigate the physical effects, we 
discuss the results of our resolution study. Numerical 
diffusivity will add to the physical diffusivity, as it is a 
natural consequence of the finite difference scheme ap- 
plied in the PLUTO code. This effect could be in par- 
ticular important in the jet launching regime close to 
the disk surface where strong gradients in density, pres- 
sure, or ma gnetic diffusivity exis t, and which may not 
be resolved. Murphy et al. (2010) have claimed that jet 
launching could be in fact present just due to numerical 
diffusivity. 

In order to estimate the impact of numerical diffusivity 
we have repeated our reference simulation case 1 with a 
three times higher resolution (case 7), but with the same 
physical diffusivity and anisotropy parameter as for the 
reference simulation. We find that the leading disk and 
outflow properties are similar to the reference case. The 
accretion rate is slightly decreased from M acc = 0.015 
(case 1) to M^c — 0.013 (case 7), while the ejection rate 
is slightly increased from M j ec = 0.008 to M c j C c = 0.009 
(see Tab. [I]). Due to the smaller computational domain 
for the high resolution simulation we cannot compare the 
asymptotic jet radii (at z — 180), however, we can com- 
pare the radius of the bulk mass flux similar to Eq. 17 fol- 
lower altitudes (at z = 60). The simulation with higher 



resolution seems to result in a slightly more collimated 
jet, with a jet radius of rj e t = 19.6 compared to fj e t = 25 
for the reference simulation. This results in a similarly 
smaller jet launching radius. The maximum jet velocities 
are just the same Uj t = 0.5. 

Figure [T4| shows the time evolution of the mass fluxes. 
We see that the accretion rates and ejection rates for case 
1 and case 7 saturate at the same level. It looks that the 
high resolution simulation needs more time to establish 
an outflow, while the accretion evolves similarly in both 
runs. One may see indication for a slightly larger accre- 
tion rate for the low resolution run, which would fit into 
the picture that the disk material can more easily diffuse 
across the field lines due to the numerical diffusivity. 

The same pi ctu re holds for the angular momentum 
fluxes (see Fig. 15), which are very similar for the kine- 
matic part, and slightly offset for the magnetic part. 



4.2. Impact of the magnetic field strength 

There is a common agreement in the literature that 
jet formation requires a certain amount of magnetic 
flux to be present in the jet launching regime. On the 
other hand, the maximum magnetic flux which can be 
supported by the accretion disk is limited by the disk 
equipartition field strength (we will neglect the question 
of the origin of the magnetic field in this paper) . 

We have also studied the impact of the magnetic field 
strength on jet launching, governed by the plasma-/? pa- 
rameter. We note that we start all simulations with the 
same magnetic field profile, however, due to diffusion and 
advection in the disk, the field distribution may change 
substantially. 

The magnetic field strength determines 

the amount of magnetic energy which is available for jet 
acceleration, and is directly interrelated with the length 
of the lever arm of the ma gnetic torque on the disk (e.g. 



Pelletier fc Pudritz| fl992| ). 

Due to the larger magnetic torque in case of a strong 
field, the disk angular momentum could be removed more 
efficiently. In order to investigate these effects quan- 
titatively, we will compare simulations with different 
plasma-/? such as case 1 with initial j3 — 10, case 8 
with f3 = 50, case 9 with /? = 250, and the weak field 
case 6 with f3 = 5000. Note that the plasma-/? is, how- 
ever, a space and time dependent function of the simu- 
lation, /? = f3(r,z,t) and not only one single parameter 
we prescribe initially at the inner disk radius. In the 
high plasma-/? regime, the flow of matter controls the 
dynamical structure of the system. For all cases of a 
weak magnetic field, the disk-jet system evolution seems 
highly perturbed - visible in the global structure of the 
outflow (Fig.pJ case 6,8,9). 

Note that the regimes with high plasma-/? are known 
to be dominated by the internal, turbulent torque, which 
is not taken into account in our simulations by an a- 
viscosity. 

In particular, for case 9 and 10, the mass flux evolution 
is exceptionally perturbed (so we decided not to include 
them in all plots and ju st s how the magnetic angular 
momentum flux). Figure 15 proves that with increasing 



plasma-/?, less angular momentum is removed from the 
disk, and, subsequently, the accretion rate is reduced. 
This is in particular visible when we compare the simu- 
lation with /? = 50 with the reference run with /? = 10. 
In case of the weakest magnetic field, thus for the simula- 
tions with /? = 250 and /? = 5000, the magnetic angular 
momentum removal is approximately negligible. For the 
weak field cases, the accretion rate decreases, and in case 
9, 10 and 6 no efficient accretion is observed. In sum- 
mary, we confirm the hypothesis that efficient magneto- 
centrifugal jet driving requires a strong magnetic flux 
(i.e. a low plasma beta), together with a large enough 
magnetic torque in order to produce a powerful jet. 

Simulation case 6 appli es a very weak magnetic field 
with /? = 5000. Figure [16] compares the two simula- 
tions case 6 with the standard setup case 7 in high res- 
olution. No smoothly structured outflow is obtained for 
case 6, but a highly turbulent outflow with non-negligible 
outflow speed and mass flux. Due to the weak poloidal 
field, the outflow in case 6 is super-Alfvenic right from 
the launching point, thus magneto-centrifugal accelera- 
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Fig. 16. — Comparison of outflow launching for low and high 
plasma-/3 cases. Display of the Alfvenic Mach number for two sim- 
ulations with three times higher resolution. The left panel shows 
the weak field case /3 = 5000 (case 6). The right panel shows the 
strong field case /3 = 10 (case 7). 

tion cannot play a role. Interestingly, the size of the 
turbulent features increases with distance from the ori- 
gin. Also the poloidal magnetic field is highly tangled 
(Fig. [l3] ) . The mass flux we measure for case 6 is about 
M e jec — 0.001 with outflow velocities of Vj ct — 0.4. While 
the velocities are comparable to case 7 (or case 1 with 
lower resolution), the ejected mass flux is substantially 
lower (factor 2). Accretion in case 6 is weak (the smallest 
of all simulations), consistent with the low angular mo- 
mentum losses (the smallest of all simulations). The ac- 
cretion rate is even smaller than the ejection rate (factor 
5), and we may call such a disk an ejection disk instead 
of an accretion disk. 

The question is what is launching and accelerating such 
a turbulent, high plasma-/? outflow? Our interpretation 
is that the initial acceleration to super-escape speed is 
by the toroidal magnetic pressure gradient (induced by 
the differential rotation between the rotating disk and 
the static corona). When the bow shock has left the do- 
main, this acceleration process decays, and the remaining 
outflow acceleration is due to weak Lorentz force (as we 
are in the super Alfvenic regime) . The vertical mass flux 
for case 6 (very weak field) is one order of magnitude 
below the mass flux in case 7 (with strong field /? = 10), 
measured at the same altitude (see Tab.flJ. In summary, 
launching conditions as in case 6 with p = 5000 do not 
produce a strong jet flow, but a light disk wind of super- 
escape speed. 

We point that the physical regime of acceleration does 
not only depend on the field strength, but also on the 
mass flux. Magnetocentrifugal effects are more evident 
for outflows with low mass flux, while in heavy jets the 
magneti c pressure gradient ma y play a substantial role 
(see e.g. Anderson et al. ( 2005[ )). 

The relative held strength (measured by the plasma- 
/?) is a function of space and evolves in time. The initial 
field and gas pressure distribution, denoted by the /?i, is 
changed by the dynamical evolution including adyection 
and diffusion of the disk magnetic field. Figure [17| shows 
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Fig. 17. — Time evolution of the plasma-/? in the inner disk-jet 
system. Snapshots of the plasma-/? distribution (in logarithmic 
scale) for easel at time = 0, 1000, 3000. 




Fig. 18. — Advection and diffusion of magnetic flux over time. 
Show is the evolution of the magnetic flux surface \& = 0.1 rooted 
initially at (r, z) = (2,0). for case 7 (high resolution) and case 
8(/3 = 50). Colors show the density distribution at t=3000 (loga- 
rithmic scale). 



the evolution of the plasma-/? for case 1. It is clearly visi- 
ble that the plasma-/? along the disk surface changes sub- 
stantially. The inner disk (for r < 15) reaches a steady 
state with a /? somewhat higher than the initial value. 
Later, around t = 3000 when the disk mass decreases 
the gas pressure decreases as well, resulting in a lower 
plasma-/? in the disk corona. The disk corona there- 
fore remains sufficiently magnetized and can support a 
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fast jet. In case 1 the magnetic field structure reaches a 
steady state balanced by field diffusion advection. 

In simulation case 7 with the same parameter setup 
but a higher resolution a steady state is reached as well. 
However, even for late time steps the magnetic flux con- 
tinues to diffuse outwards. This is somewhat surprising 
since the numerical diffusion is smaller and cannot be 
the reason (Fig. 18 top). For case 8 with the higher 
initial plasma-/?, advection seems to dominate outward 
diffusion of flux. The poloidal field distribution reaches a 
steady state with a magnetic fl ux c oncentration close to 

bottom) 



the inner edge of the disk (Fig. 18 



Thus, this 
ling as well, although 



setup seems favorable for jet launc I 
the initial plasma-/3 was high. 

As an extreme example, for simulation case 9 with the 
initial (3 = 250 the plasma-beta is considerably changed 
during the disk evolution. Case 9 shows very weak ac- 
cretion (see discussion above), so outward diffusion of 
magnetic flux dominates advection. Outward diffusion of 
flux implies a decrease of magnetic field strength. As a 
result the disk magnetization decreases, and the plasma- 
j3 reaches number values above 10 5 . Therefore this disk 
is not able to launch strong jets. So far no steady state 
is achieved. 

The interrelation between the magnetic field strength 
and o utflow launching has been discussed by other au- 
thors ( |Tzeferacos et al.||2"009{ |Murphy et al.||2010| ), in- 
dicating that regions with weak field are not able to 
generate an outflow, and that both the collimation de- 
gree and the election rate increase with stronger field. 
Murphy et al. (2010) have concluded that jet launching 
for cases of weak magnetization may be artificially sup- 
ported by numerical diffusivity within the disk surface 
layer, which should heat the gas, producing additional 
gas pressure. They suggest that ejection is possible as 
the magnetization reaches unity at the disk surface due 
to the steep density decrease. No ejection was reported 
when the mid plane magnetization becomes too small. 
Nevertheless, the asymptotic jet velocity remained too 
low to explain the observed jets speed. From our sim- 
ulations, we find that even for weak magnetization in 
the disk the disk corona is sufficiently magnetized for jet 
launching, and, depending on the efficiency of mass load- 
ing, fast jets could be driven (for comparison, see Figure 
17 for case 1). The mass loading depends on resistivity, 
and we will discuss this aspect in the next section. 

We finally consider the possibility to observe the MM 
in our simulations. This is interesting since we are deal- 
ing with magnetic fields of various strength in combina- 
tion with a differentially rotating system. To generate 
the MRI in numerical simulations, two conditions are 
essential - a large enough (but not too large) plasma- 
/3 in the disk as to initiate instability, and also a grid 
of high enough resolution in order to avoid damping of 
small wavelengths. The wavelength of most unstable 
MRI mode is given by Amri = 27rt>A,z/^K for a Ke p- 



lerian rotation J7k (see e.g. (Romanova et al. 2011)) 
When we calculate this wavelength for a simulation with 
/3 = 10, we obtain Amri — 0.057, implying that we can- 
not resolve the MRI with our setup as the grid size is 
about half of the MRI wavelength. 

We point out that the magnetic diffusivity applied in 
our simulations also contributes to suppress the genera- 
tion of MRI turbulence. 
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Fig. 19. — Mass flux M c j oc vs poloidal diffusivity. Shown is the 
ejected jet mass flux (triangles) and the accretion rate (circles) as 
function of the poloidal diffusivity for simulations with the same 
magnetic field strength, i.e. case 1 (r; Pi i = 0.03), case 2 (r; Pj i = 
0.01), case 3 (r? Pii = 0.15), case 4 (r? Pji = 0.09), case 5 (7? Pii = 0.03). 



4.3. Magnetic diffusivity strength and anisotropy 

The strength and anisotropy of the magnetic diffusiv- 
ity - controlled by the coefficients rj^^r/pi and \ ( see 
Tab. [T]) - are essential parameters governing the coupling 
between the magnetic field and the plasma. Efficient 
magneto-centrifugal driving requires strong coupling be- 
tween the magnetic field and the rotating disk material, 
thus a rather low diffusivity. A similar argument holds 
for the launching by magnetic pressure gradient. On the 
other hand, mass loading from accretion to ejection re- 
quires a certain degree of diffusivity in order to transport 
mass across the field lines. 

So far, no general model for the magnetic diffusivity 
distribution in disk-jet structures is available. In many 
simulations considering magneto-centrifugal driven out- 
flows, the outflow is governed by ideal MHD, while the 
diffusivity is confined to the disk. In order to investigate 
the impact of diffusivity on jet launching, we have there- 
fore performed simulations with varying strength of the 
diffusivity components r\^^ and 7y Pj ;. 

We first compare simulations with the same anisotropy 
parameter \ = 3 as reference run case 1. In general, the 
strength of magnetic diffusivity governs the disk accre- 
tion rate. Our reference run case 1 reaches quasi steady 
state, establishing a balance between inward advection 
and outward diffusion of magnetic flux. This happens 
after t = 1000 and is disturbed again in the late stages 
of the dynamical evolution, most probably due to the 
change in the overall disk dynamics due to the decreas- 
ing disk mass. A change in diffusivity will also affect this 
balance. We now compare case 4 with higher and case 
2 with lower diffusivity. In the less diffusive case, ad- 
vection dominates and we obtain a higher accretion rate. 
We find that not only the accretion rate, but also the 
ejection rates are higher for disks with lower diffusivity. 
This trend is shown in Fig. 
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where we display the ac- 
cretion rate and also the ejection rate as function of the 
the poloidal diffusivity. 

Two physical processes affect the disk dynamics - ad- 
vection and diffusion. In a less diffusive disk, the plasma 
is stronger coupled to the field, and advection of magnetic 
flux dominates. Consequently, the magnetic flux surfaces 
are located further in the inner part of the disk, "ro- 



18 



Shcikhnczami ct al. 



1.0 



0.3 



r, 



4 



0.2 



0.0 1 

0.00 



0.05 



0.10 



0.15 



0.20 



Fig. 20. — Jet velocity vs poloidal diffusivity at t = 3000. Shown 
is the typical jet velocity as a function of but different diffusivity, 
i.e. case 1 (r] p ^ = 0.03), case 2 {rj p ^ = O.Of ), case 3 (r; p i = 0.15), 
case 4 (»7 P; j = 0.09), case 5 (r? Pi i = 0.03). 

tate" faster, and a stronger is induced. The stronger 
toroidal field may lead to a stronger acceleration of the 
outflow, by stronger Lorentz forces (either magnetic ten- 
sion of pressure forces). In addition, by comparing the 
vertical profiles of net launching forces at t — 1000 for 
less diffusive disk (case 2), and reference run (case 1), we 
find a larger net force in case 2 by which implies that 
more disk material can be lifted into the outflow. 

This is exacly what we observe in the mass fluxes - 
the ejected mass flux in case 2 is two times the ejected 
mass flux in case 1. Note that a stronger will impose 
stronger pinching forces ~ (B$ ■ V)B^ on the disk and 
the ejected material, thus opposing the loading process. 
However, the magnetic pressure gradient of the toroidal 
field ~ V-B^ has a radially outward direction, thus sup- 
porting ejection. By comparing the vertical profiles of 
both terms for the two cases with f3 = 10 and f3 = 50, 
we find that i) the toroidal magnetic field pressure gra- 
dient component in vertical direction (VBV) Z is always 
dominating the pinching force of the toroidal field, and 
that ii) for higher plasma-/3 the (V£?^) 2 is larger than 
for the lower plasma-/3 case. Thus, in our simulations 
the toroidal field is supporting ejection. 

In summary, although a certain magnitude of diffu- 
sivity is required for mass loading, the lower diffusivity 
allows for enhanced mass loading of accreting material 
into the outflow. Consequently, less diffusive disks tend 
to have higher ejection rates . This result is in a greement 
with the previous literature 



Zanni et al. (2007) 



Our simulations confirm that the strength of magnetic 
diffusivity does affect the asymptotic speed of the out- 
flow. Figure [20| shows the typical jet velocity versus the 
poloidal diffusivity. Comparing different simulation runs 
(see Tab.PLJ), we find that the outflow from the less dif- 
fusive disK (case 2) remains slower, while the outflows 
formed from more diffusive disks (case 3 and 4) are accel- 
erated to higher speed. For case 2 with three times lower 
magnetic diffusivity the asymptotic speed is reduced by 
ps 30%, while for case 4 with three times higher magnetic 
diffusivity the increase is about ps 30% (see Tab. [I). This 
is reasonable from a physical point of view, since in a 
weakly diffusive disk more material is diverted into the 
outflow. Thus, with the same amount of magnetic flux 
available, correspondingly less magnetic energy could be 
transferred per outflow mass unit. The more massive 



outflows are only accelerated to lower speed if launched 
with similar magnetization. Alternatively, a higher mag- 
netic diffusivity resulting in a weaker mass load, leads to 
a faster outflow. 

When comparing our results to the previous literature, 
one should keep in mind that we have applied a mass flux 
weigthed velocity which we consider as the typical veloc- 
ity of the bulk mass flux. The velocity value is generally 
lower as the maximum speed we measure in the outflow 
and which is mostly given in the literature. While the 
typical speed ranges from about 0.5 to 0.9 (in code units), 
the maximum speed in the outflow ranges from 1.2 to 
1.8 inner disk Keplerian velocities and is comparable to 
other paper in the literature. Furthermore, we have de- 
tected a slight time evolution in the velocity (see also the 
mass flux evolution shown in Fig. 
run easel, 



14). In our reference 
the maximum asymptotic speed varies from 



u p,max 

comparison, 



1.8 at t = 500 to v n 



Zanni et al. (2007); Tzeferacos et al 
give a maximum speed of f p , m ax 



= 1.2 for t > 1000. For 



(2009) 



1.5.. .4.5 at t 

Next, we investigate a possible impact of anisotropy 
in the diffusivity-tensor, parametrized by \. We com- 
pare case 1 and case 5 with the same polodial diffusivity 
r/ p = 0.03, but with case 5 having a lower diffusivity in 
toroidal direction. By comparing their leading properties 
(see Tab. Ill, we find them quite similar properties. In 
particularTfor both cases (with the same poloidal diffu- 
sivity) the ejection rates are also similar. However, case 5 
with less toroidal diffusivity has a higher accretion rate. 

Anisotropy of magnetic diffusivity does also impact the 
rotation of the outflow. It is believed that jets are rotat- 
ing and the rotation is basically driven by the underlying 
disk rotation and the coupling between matter and field 



dBacciotti et al.|2002[ [Anderson et al.|2003 
20041 IFcndt 201 lk. For a high toroidal d 



Coffey et al 
mifusivity, the 
coupling is weak and, thus the acceleration in toroidal 
direction is weak as well. With our simulations we can 
approve this concept. Figure 21 shows the rotational ve- 
locity for the lower part of the jet. The outflow resulting 
from simulation case 1 with a nine times higher toroidal 
diffusivity shows a 50% lower rotation rate than for case 
5. Similarly, the outflows in simulation case 4 with a 
three times higher toroidal diffusivity shows a 50% lower 
rotation rate than case £ p"^l 

We observe a close correlation between the, accretion 
rate M acc and angular momentum flux (Jkm + ^mag) from 
the disk. In a system with higher angular mom entu m re- 
moval, a higher accretion rate is observed (Fig. 15 ), with 
the accretion rate in case 5 being larger than tor case 
1, and for case 1 larger than for case 3. This confirms 
the common believe that in order to obtain higher accre- 
tion rates, a more efficient angular momentum removal 
is required. 

To discuss the interrelation between collimation and 
diffusivity, we again apply a collimation degree by com- 
paring the directed mass fluxes (see Eq. 18 ). As a general 
result we find that outflows launched from disks of higher 
poloidal diffusivity rj v are more collimatcd. This interre- 
lation are displayed in Fig. [22j 

We see case 2 with the magnetic diffusivity of r) p ; = 

13 We cannot compare cases 1 and 4 or 5 and 4 directly, as, due 
to their different poloidal diffusivity these outflows have different 
mass fluxes 
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Fig. 21. — Outflow rotation and toroidal magnetic diffusivity. Shown is the rotational velocity vj, of the inner disk-jet system at t= 3000 
for simulation runs with different toroidal magnetic diffusivity, respectively different anisotropy parameters \, i.e. cases 1, 5, 3, 4 (from left 
to right). 
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Fig. 22. — Collimation degree and jet radius vs poloidal diffusiv- 
ity. Shown is the collimation degree £ (filled triangles) and the mass 
flux weighted jet radius rj et (filled circles) over the poloidal diffusiv- 
ity for simulations with the same magnetic filed strength, but and 
different diffusivity, i.e. case 1 (%>.i = 0.03), case 2 (r) Pt i = 0.01), 
case 3 {rjp,i = 0.15), case 4 (r) Pj i = 0.09), case 5 (?7 Pi i = 0.03). 

0.01 - maybe below a critical value, evolves differently 
from the others and has both a very low degree of col- 
limation and an exceptionally large mass flux weighted 
jet radius. The inner part of this outflow appears rather 
hollow (Fig. [13]). 

The comparison of the (mass flux weighted) jet radii 
for these runs provide a similar picture. We find an 
interrelation such that for increasing magnetic diffusiv- 
ity the jet radius decreases such that rj e t(case 3) < 
r jet ' 



(case 4) < rj Ct (case 5) < rj ct (case 1) < rj et (case 2) 



4.4. Scale height of magnetic diffusivity 

We have argued above that the disk material which 
is launched to form an outflow is likely to be turbu- 
lent, so we may expect a magnetically diffusive disk wind 
above the disk surface. Thus, since mass loading requires 
poloidal magnetic diffusivity, the jet launching area can 
be extended into higher altitudes above the disk surface. 
Furthermore, the jet launching area is extended into a 
domain where the plasma-/? is comparatively low. Tak- 
ing into account both effects, one may therefore expect 
i) to launch jets from disks which would themselves be 
insufficiently magnetized for jet-driving, and also ii) to 
launch outflow which have a rather low mass flux and 
subsequently higher terminal velocity. In order to inves- 
tigate this effects, we have therefore prescribed different 
height scales H, n fo r the magnetic diffusivity in our simu- 
lations (see Eq 13 1 . The physical properties measured for 



14 For simulations with higher resolution (case 6 and case 7), we 
adopt a smaller physical grid size with r ou t = 50, so we cannot 
compare the jet radius properly. 



TABLE 2 

Comparison between different diffusive scale heights. 
Shown are some of the physical properties mentioned in 
Tab. [l] for different diffusive scale height runs which 
have been carried out up to t=2000. 





case 1 


case 11 


case 12 


case 13 




0.1 


0.2 


0.3 


0.4 


Ar 


0.064 


0.064 


0.064 


0.064 


As 


0.066 


0.066 


0.066 


0.066 




0.03 


0.03 


0.03 


0.03 


n<t,,i 


0.09 


0.09 


0.09 


0.09 


X 


3 


3 


3 


3 


P 


10 


10 


10 


10 


r jct,z=180 


43.8 


73.3 


56.4 


48.4 


Jjet,z=60 


25.01 


45.8 


39.2 


34.4 


n 


3.8 


7.2 


7.2 


6.2 


"jet,z=170 


0.46 


0.35 


0.33 


0.41 


l>jet,z=280 


0.53 


0.46 


0.41 


0.47 




0.015 


0.024 


0.022 


0.018 


M=jcc 


0.005 


0.015 


0.015 


0.01 


^kin 


0.01 


0.005 


0.009 


0.01 




0.033 


0.015 ±0.01 


0.02 ±0.003 


0.028 



different diffusive scale heights are shown in Tab. [2] In 
general, we find different stages in the temporal evolution 
of the mass flux (accretion and ejection) in the disk-jet 



system (Fig. 23). For small diffusive scale heights H^, a 
turbulent pattern and perturbations appear in the out- 
flow resulting in a more filamentary outflow. It seems 
that these perturbation have larger amplitudes in case of 
a smaller diffusive scale height. The amplitudes decrease 
when the diffusive scale height decreases. For e„ = 0.4 
they are damped completely. We believe that these per- 
turbations are generated by a physical process within the 
launching region and that they are simply smoothed out 
by the magnetic diffusivity. A further study is needed to 
clarify this issue. 

For increasing scale height the accretion rate and ejec- 
tion rate decrease, confirming the correlation derived 
above between the increasing magnitude of diffusivity 
and decreasing mass fluxes. In other words, an increasing 
diffusivity scale height simply corresponds to a higher dif- 
fusivity with the consequences as discussed above. Inter- 
estingly, all our simulations for different diffusivity scale 
heights converge to the same kinetic angular momentum 
flux in the outflow. On the other hand, we find a trend 
of decreasing magnetic angular momentum flux with de- 
creasing diffusivity scale height. 



20 



Shcikhnczami ct al. 




o.os 



o 01 



0.00 
2000 



0.04 - 




Fig. 23. — Effect of the magnetic diffusivity scale height. Shown is the evolution of the accretion and ejection rate Mkim Af mag (top), 
and the evolution of the kinetic and magnetic angular momentum flux J^in, -/mag (bottom), for simulations applying a different diffusive 
scale height H v = 0.1 (solid line), 0.2 (dotted line), 0.3 (dashed line), 0.4 (dot-dashed line). 



along with the evolving outflow, the magnetization above 
the disk surface is substantially lower compared to the 
initial condition. 

We obtain jet velocities which seem to be smalle r com- 
pared to pre viously published values. For example, |Zanni| 



5. COMPARISON WITH PREVIOUS NUMERICAL STUDIES over 1000s of dynamical time scales (Fig. |i7|). Thus 

In this section we highlight some of our results dis- 
cussing them in the context of numerical studies pub- 
lished previously. Our main goal was to perform a consis- 
tent parameter study to investigate a variety of physical 
effects involved in jet launching by using one code with 
one numerical setup. Only this allows us to disentangle 
the leading effects involved in jet launching without the 
uncertainty introduced by interpreting results obtained 
with different codes (VAC, PLUTO, FLASH) or setups. 

We first stress the point that we apply a well defined, 
fixed-in-time magnetic diffusivity distribution. The mag- 
netic diffusivity is an leading parameter which infl uences 
the whole disk-jet evolution flZanni et a l. 2007). Pre- 
vious simulations ( Casse fc Keppcns 2002 



et al. (20071 detect asymptotic velocities of 1.5-1.8 times 
the footpomt Keplerian speed for the field lines crossing 
the upper boundary. Simulations of jet forma tion from 



a fixed-in-time disk surface boundary (see e.g. |Ouyed fc 



Pudritz 
values 



Zarm i et al. 

2007| |Tzeferacos et al.||2009[ |Murphy et al.||2010[ ) have 



applied a ditlusivity evolving in time. This might be 
considered as more advanced, however, the previous au- 
thors do not provide details about time-evolution of the 
diffusivity, thus intruducing some level of uncertainty. 
Here, we have demonstrated how diffusion and advec- 
tion of magnetic flux does change the disk Alfven speed 
(and thus the magnetic diffusivity if a time-dependent 
77 (t) ~ VA(t) is applied) by a factor of ten. On the other 
hand, by applying different strength and scale height for 
the diffusivity we investigate the interrelation between 
magnetic diffusivity and the accretion and ejection mass 
fluxes. 

Similarly, we provide information about the spatio- 
temporal evolution of the magnet ization (or plasma-/3 in 



our case). |Tzeferacos et al. (2009J have published a mag- 
netization stlicIyTTioweverTon a smaller grid and for sub- 
stantially shorter dynamical evolution. In our paper we 



( |1997| ); |Fendt fc Cemeljic| ( ]2002| ) ) give similar 
In contrast, we report lower velocities of about 
0.8 times the Keplerian speed at the inner disk radius. 
Note however, that our values result from a differently 
defined measure of the jet velocity. We suggest that this 
mass-flux weigthed asymptotic speed - the average veloc- 
ity of the bulk mass flux - is more applicable to obser- 
vations. The maximum jet velocities we observe in our 
simulations are higher and also up to 1.8 time the Kep- 
lerian speed at the inner disk radius, which is consistent 
with the previous studies. We explicitely consider the jet 
rotation in our simulations a topic which has not been 
treated in previous launching simulations. 

We quantify the jet collimation degree by the directed 
mass fluxes. For the jet opening angle we find somewhat 
larger values than in previous studies. We suppose that 
this results from the modified outflow condition we have 
applied as we have seen that the original outflow condi- 
tion lead to smaller opening angles. 

We have run long-lasting simulations of up to 5000 dy- 
namical time-scales U — ri/vK,i- This allows us to i) 
reach a quasi steady-state situation of the simulation, 
but also to investigate the evolution beyond. We find a 



show in particular how the magnetization profile flares sli S nt but persistent change in the inflow-outflow dynam- 
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ics on this very long time scales, such that the accretion 
rate slightly decreases. The main reason for this we see 
in the decrease of the disk mass for such long time scales. 
However, the effect has an direct astrophysical applica- 
tion as it may be applied to the long-term evolution of 
Classical TTS when the disk accretion in fact weakens 
after some 10 6 years. We note that previous simulations 
were substantially shorte r, lasting until 400 dynamical 

]T zeferacos| 
Kcppens 



2004 
ByP 
isk 



time steps ti — n/wK.i ( Zanni et aLN 2007 
et al.|2009) or 30 inner d isk rotations (Uasse 
200G 



Meliani et al. 



An exception is the simulation 
urphy et al. 1 2010p lasting for more than 900 inner 
disk orbits. That paper focusses, however, on the launch- 
ing physics and does not show the large-scale evolution 
of the outflow. 

Together with the long time evolution, w e have also 
have applied a large grid size. Compared to |Tzeferacos| 



et al.] (2009[ ) who investigate how magnetization affects 
launching by using a 40 x 120 n grid , our grid has about 
double size. Murphy et al. ( 2010[ ) concentrate on the 
launching process of weakly magnetized disk and display 
only results on 40 x 40 n images (which are a subsets of 
a 280 x 840 ri). For YSO, our grid extension corresponds 
to about 28 AU along the jet, well into the observable re- 
gion. We note that we have applied an equidistant grid 
for the magnetically diffusive disk area which is consis- 
tent with the PLUTO code requirements. This is differ- 
ent to Murphy et al. ( |2010 ) who attach a scaled grid for 
the diffusive outer disk! 

6. CONCLUSIONS 

We have presented results of MHD simulations investi- 
gating the launching of jets and outflows from a magnet- 
ically diffusive disk in Keplerian rotation. The time evo- 
lution of the accretion disk structure is self-consistently 
taken into account. The simulations are performed in ax- 
isymmetry applying the MHD code PLUTO. The main 
goal of our simulations was to study how magnetic dif- 
fusivity (its magnitude and distribution) and magnetiza- 
tion affect the disk and outflow properties, such as mass 
and angular momentum fluxes, jet collimation, or jet ra- 
dius. 

Our grids extend to (96 x 288) inner disk radii with 
a resolution of (0.064 x 0.066), respectively (50 x 180) 
inner disk radii with a higher resolution of (0.025x0.025). 
An internal boundary (sink) is placed close to the origin 
absorbing the accreted mass and angular momentum. 

We have prescribed a magnetic diffusivity in the disk 
based on an a-prescription. One of our parameters was 
the scale height of the magnetic diffusivity with the op- 
tion to have it higher than the thermal scale height. This 
can be motivated considering that it is the turbulent disk 
material which is loaded into the outflow, and that the 
turbulence pattern is swept along with the disk wind un- 
til it decays. We have investigated disks carrying a mag- 
netic flux corresponding to an initial plasma-beta ranging 
from 10 to 5000 at the inner disk radius. 

As a general result we observe a continuous and robust 
outflow launched from the inner part of the disk, expand- 
ing into a collimated jet and is accelerated to super fast 
magneto-sonic speed. The key results of our simulations 
can be summarized as follows. 

(1) Concerning the acceleration of the outflows, our 
simulations confirm that the magneto-centrifugal accel- 



eration mechanism is most efficient in the low plasma-/? 
regime, while for weak magnetic fields the toroidal mag- 
netic pressure gradient drives the ejected material. 

We also confirm that the magnetocentrifugal mecha- 
nism also depends on the mass load in the outflow, as 
this mechanism works more efficiently for outflows with 
low mass fluxes. However, compared to the magnetic 
pressure driven outflows, jets in the magnetocentrifugal 
acceleration regime have usually higher mass fluxes. 

In our simulations with very high plasma-/3 we detect 
a highly unsteady behaviour 

(2) Efficient magneto-centrifugal driving which can ac- 
celerate jets to high kinetic energy relies on a strong cou- 
pling between magnetic field and the rotating disk, thus 
on a low diffusivity. We find that it is the poloidal diffu- 
sivity ?7 P which mainly affects the driving of the outflow. 
However, besides the coupling needed for acceleration, 
also the launching of material depends on diffusivity. 
With increasing rj p , the mass fluxes (both the accretion 
rate and the ejection rate) decrease. Subsequently, the 
higher ejection rates result in a lower asymptotic outflow 
velocities. 

(3) We measure typical outflow velocities are in the 
range of 0.3 — 0.8 times the inner disk rotational ve- 
locity with the tendency that the mass fluxes obtained 
in magneto-centrifugally driven outflow are substantially 
higher. Here, we confirm the clear (inverse) correlation 
between jet velocity and mass load, as it is well known 
from the literature. Note that we have apply a mass- 
flux-weighted jet velocity which we find more applicable 
to the observations. For the bulk mass flux we find lower 
velocities compared to other papers, which are mostly 
dealing with the maximum speed obtained in the simula- 
tion. Our maximum velocities are similar. The relatively 
high speed of the outflows with low Poynting flux which 
are driven by poloidal pressure gradient is unclear. 

We find that the toroidal diffusivity affects the outflow 
rotation - a small toroidal diffusivity implies a larger jet 
rotational velocity. This has not shown before in simu- 
lations. 

(4) We do not find a clear correlation between the out- 
flow collimation and the magnetic field strength. Also 
weakly magnetized outflows, which are driven by the 
magnetic pressure gradient, and which we find to be quite 
unsteady, show a high degree of collimation. The ques- 
tion of collimation for the weakly magnetized outflows is 
not answered. 

We find that outflows within the magneto-centrifugal- 
driving regime the flows ejected from a weakly diffusive 
disks are only weakly collimated. Similarly, their jet ra- 
dius (here defined as mass flux weighted radius) is larger 
in case of a lower poloidal magnetic diffusivity. 

Follwoing the magnetic flux surface along the bulk 
mass flux from the asymtotic regime to the launching 
area, we can defined the launching area of the outflow. 
We find a size of the launching area from which the bulk 
of the mass flux originates in the range between 3 and 8 
inner disk radii or about 0.4 AU. 

(5) Depending on the strength of magnetic diffusivity, 
the disk-jet structure may evolve into a steady-state. We 
found that the cases with the strong field with $ ~ 10 
and poloidal diffusivity r/ p ^ > 0.03 will reach a quasi 
steady state, confirming the literature. 

(6) The magnetic flux profile along the disk is subject 
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to advection and diffusion. We find that the magneti- 
zation (or plasma-/?) of disk and outflow may therefore 
substantially change during the time evolution. We have 
observed that the initial disk magnetization may change 
by a factor of 100. This may have severe impact on the 
launching process and the formation of the outflow in the 
sense that a rather highly magnetized disk may evolve 
into a weakly magnetized disk which cannot drive strong 
outflows. This issue has not been discussed before in the 
literature. 

(7) For very long time scales the accretion disk changes 
its internal dynamics, as due to outflow ejection and disk 
accretion the disk mass decreases. As a consequence, the 
accretion and ejection rates slightly decrease. In order to 
compensate for this effect, we have applied a large outer 
disk radius providing a large mass reservoir for the inner 
jet-launching disk. 

(8) For our simulations, we found that 10-50% of the 
accreting plasma can be diverted into the outflow. For 
i) less diffusive disks, ii) a strong magnetic field, iii) a 
low poloidal diffusivity, or a iv) lower numerical diffu- 
sivity (resolution), the mass loading into the outflow is 
increased - resulting in more massive jets. We inter- 
pret as physical reason the more efficient extraction of 
angular momentum from the disk, due to the stronger 
matter-field coupling. Note that we do not consider in 
our simulations viscosity or the wealth of t hermal effects 



which play an essential role for launching ( Casse & Fer- 
"12000 i 



reira 



(9) We found that jets launched in a setup with smaller 
diffusive scale height are more perturbed. The same ef- 
fect is seen in outflows launched from disks with weaker 
diffusivity. 

(10) We finally remark that it is essential to do long- 



term simulations covering thousands of rotational pe- 
riods in order to find a steady state situation of the 
accretion-ejection dynamics. Our simulations run for 
5000 dynamical time steps, corresponding to about 900 
revolutions at the inner disk radius as adopted in our 
reference run. 

In summary, we confirm the hypothesis that efficient 
magneto-centrifugal jet driving requires a strong mag- 
netic flux (i.e. a low plasma beta), together with a large 
enough magnetic torque in order to produce a powerful 
jet. In addition, the magnitude of (turbulent) magnetic 
diffusivity plays the major role in the ejection efficiency, 
while the anisotropy in the diffusivity mainly affects the 
jet rotation. Both results imply that the structure of the 
asymptotic jet is indeed governed by the properties of the 
accretion disk, here parametrized by the magnetization 
and magnetic diffusivity. The mass ejection-to-accretion 
ratio along with the momentum and energy transfer rates 
from inflow to outflow are essential properties for any 
feedback mechanism in star formation or galaxy forma- 
tion scenarios and could only be derived from simulations 
resolving the inner region of the jet-launching accretion 
disk. 
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APPENDIX 

UNITS AND NORMALIZATION 

We normalize all variables, namely r,p,v,B to their fiducial values at the inner disk radius. We have adopted the 
following typical number values for a YSO of M — 1 M and an AGN of M = 10 s M . A change in one of the system 
parameters accordingly changes the scaling of our simulation results which can thus be applied for a variety of jet 
sources. The inner disk radius is usually assumed to be a few radii of the central object, 

r, = 0.028 AU ( ( (YSO) 

(AGN), (Al) 



3-Ryso/ \ / 



M 



V 10i? s / V 10 8 M G 



where R$ = 2GM/c 2 is the Schwarzschild radius of the central black hole. The inner disk radius is usually assumed 
to be located at the marginally stable orbit at 3 Rg ■ Since we apply the non-relativistic version of the PLUTO code, 
we cannot treat any relativistic effects. We therefore apply a scaling of n — lORg. For simplicity and convienient 
comparison with the previous literature we apply r\ — 0.1 AU for most comparisons concerning stellar sources. The 
orbital velocity at the inner disk radius is 

- 6 - 7xl ° 4kmS "(rot) (lO^o) (AGN). (A2) 

For a r; = 0.1 AU distance from a YSO the orbital speed is Uk,i = 94kms~ 1 . The mass accretion rate is a parameter 
which is in principle accessible by observation. Subject to the disk model applied, the observed disk luminosity can 
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be related to an accretion rate. For a YSO, the accretion rate is typically of the order of M acc ~ 10 _7 M Q yr~ 1 , 
providing a normalisation of the density p\ with M\ = rfpiVK.i Applying a length scale r; = 0.1 AU and a velocity scale 
^k.i = 94kms , we obtain 



M i = lO- 5 M yr- 1 
= lOM yr" 1 



M 
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Normalized value of the magnetic field is obtained by considering the plasma-/? and the field strength at the equator 
at the inner radius By = y / 87rH7A, 
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A MODIFIED RADIAL OUTFLOW BOUNDARY CONDITION 



An essential point of our setup is to impose an outflow boundary at the outer part of the domain in r-direction 
avoiding artificial collimation. We have implemented a current-free outflow boundary condit ion which prevents spurious 
collimation by Lorentz for ces and which has been thoroughly tested in our previous papers ( Porth fc Fendt|2010 Porth 
et al.|[20lT Vaidya et al.||20lT I . The boundary condition considers a vanishing toroidal electric current density in the 
ghost cells. Therefore, the Lorentz force component perpendicular to the grid boundary vanishes. The boundary 
constraints are applied on the poloidal and toroidal component of the electric current density, j r — —d z B<p, j z — 
r~ 1 d r rB ( p, and = d z B r — d r B z , respectively. 

In principl e, this is considered f or the grid ghost cells (i C nd, j), adjunct to the domain boundary at (i cn d + 1, j) (see 
the §3.1.1 in ( |Porth fc Fendt|2010[ )). The transverse B t (i end ,j + 1/2) and normal B n (i cnd + 1/2, j), B n (i cnd + 1/2, j + 1) 
magnetic field components in the domain, together with the transverse field component B t (i cn d + 1, j ' + 1/2) of the 
first ghost zone, constitute a toroidal corner-centered electric current I^,(i en d + 1/2, j + 1/2). With that a constraint 
for a current-free 1^ = boundary condition can been implemented numerically 
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